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We develop a non-equilibrium theory to describe weak Coulomb blockade effects in open quan- 
tum dots. Working within the bosonized description of electrons in the point contacts, we expose 
deficiencies in earlier applications of this method, and address them using a 1 /N expansion in the 
inverse number of channels. At leading order this yields the self-consistent potential for the charging 
interaction. Coulomb blockade effects arise as quantum corrections to transport at the next order. 
Our approach unifies the phase functional and bosonization approaches to the problem, as well as 
providing a simple picture for the conductance corrections in terms of renormalization of the dot's 
elastic scattering matrix, which is obtained also by elementary perturbation theory. For the case of 
ideal contacts, a symmetry argument immediately allows us to conclude that interactions give no 
signature in the averaged conductance. Non-equilibrium applications to the pumped current in a 
quantum pump are worked out in detail. 

PACS numbers: 73.23.-b,73.21.La,73.23.Hk 



o 

> 
00 

o 
in 
o 



I 

o 
o 



I. INTRODUCTION 

Coulomb blockade is the phenomenon that transport 
through an almost isolated system is prohibited by the 
energy cost Ec to add or remove an electron. The 
Coulomb blockade can be lifted by fine tuning a gate 
voltage to a point of charge degeneracy, where the en- 
ergy cost of adding one electron is zero, or by raising 
the bias voltage or the temperature above Ec- A fun- 
damental question is to what extent the quantization of 
charge on the system, and hence the Coulomb blockade, 
is affected by the inclusion of that system in an electrical 
circuit. This question has been addressed both for metal 
particles coupled to electrodes via tunnel barriers, as a 
function of the dimensionless conductance of the tunnel 
barrier, and for semiconductor quantum dots coupled to 
electrodes via ballistic point contacts, as a function of the 
number of propagating channels in the point contacts. In 
each case, the Coulomb blockade is lifted when the con- 
ductance of the contact to the electrodes is larger than 
the conductance quantum /h. Such quantum dots (or 
metal particles) are referred to as 'open', in contrast to 
'closed' quantum dots, which are connected to the elec- 
tron reservoirs through tunneling contacts with a con- 
ductance smaller than the conductance quantum. The 
present paper deals with Coulomb blockade in open semi- 
conductor quantum dots, which we shall refer to as 'weak 
Coulomb blockade'. 

In the literature, two approaches have been taken to 
this problem. One theory of 'weak Coulomb blockade' 
was proposed by Flensberg^ , MatveeviS, and Furusaki and 
Matveevf2i4 who realized that a tractable description of 
the strong quantum charge fluctuations can be formu- 
lated using the one-dimensional nature of the point con- 



tacts. Whereas these theories accounted for the electron- 
electron interactions non-perturbatively, they neglected 
the effect of electrons coherently traversing the dot or 
being backscattered from inside. Extensions of the the- 
ories of Refs. [l|Qy|j to the case of open quantum dots 
with coherent scattering of electrons were given by Yi 
and Kane^ for quantum dots in the quantum Hall regime 
and by Aleiner and Glazman^ and Brouwer and Aleiner— 
for the general case (see Ref. |3 for a review). 

The other theoretical approach starts from an effective 
theory in which the primary dynamic variables are the 
potential differences between the quantum dot and the 
electrodesi^ This approach has mainly been used to study 
metal grains coupled to the reservoirs via tunnel bar- 
riers with many Channels,i"i"ii2,i3,i4,i5,i6,i7,i8,i9,20^ ^ut 

recently it has been applied to semiconductor quantum 
dots, both withoutSiiSSiS^ and with24 coherent scattering 
of electrons from inside the quantum dot. The general 
conclusions of both approaches are the same: For the in- 
coherent case and ideal point contacts (no backscattering 
in the contact), it was found that all charge quantization 
effects completely vanishes iSiiiiSi In all other cases — 
contacts with a small amplitude Tc for backscattering in 
the contact or coherent scattering from inside the dot — 
, charge quantization gives a correction to the dot's ca- 
pacitance or conductance. This correction is small, but 
it becomes more important as the temperature is low- 
ered or the number of channels in the point contacts is 
reduced. 

The previous works on 'weak Coulomb blockade' all 
dealt with thermodynamic properties (capacitance) or 
time-independent transport. It is the goal of the present 
work to develop the general non-equilibrium theory for 
Coulomb blockade in open quantum dots, accounting for 
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coherent scattering of electrons inside the dot. Our main 
result is an expansion in powers oi\/N for the current 
through the dot, where N is the total number of channels 
in the point contacts. 

Two examples of non-equilibrium problems involving 
quantum interference in open and interacting mesoscopic 
systems have been a particular motivation for this work. 
First, the possibility of applying time dependent poten- 
tials via local gating of quantum dots led to the theoret- 
ical and experimental investigation of 'quantum pumps', 
in which a periodic perturbation of potentials inside the 
quantum dot, combined with quantum interference, leads 
to a dc current through the dot»2S, For quantum dots 
with time dependent potentials, effects of coherent scat- 
tering inside the quantum dot, charge quantization, and 
electron-electron interactions were considered s eparately 
theoretically (see, e.g., Refs. I26l2i28l29l30l3l|) . but not 
together. The second example that motivated this work 
is the experiment of Ref. 132 . where the nonequilibrium 
electron distribution function of a current carrying diffu- 
sive wire was measured by tunneling spectroscopy. This 
allowed an experimental determination of the electron- 
electron collision integral, and led to the discovery of a 
new mechanism of energy relaxation 

Following previous studies of Coulomb blockade in 
open quantum dots, the only electron-electron interac- 
tion term we consider is the capacitive charging energy 
of the quantum dot. This requires relatively large quan- 
tum dots: For open quantum dots Coulomb blockade ef- 
fects are small in the inverse of the number of channels N 
in the contacts connecting the dot to the electron reser- 
voirs, whereas the residual parts of the Coulomb interac- 
tion are small in the inverse dimensionless conductance 
of the 'closed' dot&^iS^ The ratio of these two conduc- 
tances equals the ratio of the dot's dwell time and ergodic 
time, which is large for large quantum dots. This param- 
eter justifies the use of the random matrix description 
of scattering by the dot, and means that the problem is 
effectively 'zero dimensional'. The dimensionality of the 
system appears only in the corrections to this picture. 

Unlike for the case of linear time-independent trans- 
port, where the capacitive interaction gives corrections 
to transport properties beyond the Hartree level only, 
interaction corrections to nonlinear or time-dependent 
transport exist already when interactions are described 
by a self-consistent (Hartree) potential Vd. As pointed 
out by Biittiker and coworkers, the reason is that, in 
non-equilibrium or time-dependent problems, a change 
in bias voltage or a change in a gate voltage potentials 
may change the number of electrons on the dot, which, 
in turn, changes VdrS*2L2& 

^d(0 = -^((^dot(i))-AA). (1.1) 

Here C is the dot's capacitance, — e(iVdot) is the dot 
charge, and —eM the offset charge induced by the gates. 
The effect of a time-dependent or bias-dependent Hartree 
potential Vd on mesoscopic fluctuations of various trans- 



port properties — but neglecting the 'weak Coulomb 
blockade' — was investigated in Refs. I30il36ll3'i1138ll39l 
To leading order in 1/iV, our theory recovers the self- 
consistent theory with the Hartree potential Hl.l|l . In 
this sense, our result can be viewed as a formal confir- 
mation of the self-consistent theory. Further, it gives 
a controlled method to find interaction corrections that 
cannot be described by means of a self-consistent poten- 
tial. 

A second motivation of this work is to compare and 
unify the results from two different approaches to these 
corrections — the 'one-dimensional' and 'environmental' 
formalisms. Whereas both approaches have given equal 
results for the case of quantum dots without coherent 
scattering from inside the dot, results for the dc con- 
ductance of coherent dots, as reported by Brouwer and 
Aleiner— and Golubev and ZaikinS4, are in disagreement. 
Using a formalism closely related to that of Ref. 0, we 
are able to present detailed calculations in two limits: 
A formal expansion in the scattering matrix of the dot 
for an arbitrary number of channels TV in the contacts, 
and as an expansion in powers of 1/N . The former limit 
reproduces the result of Ref. 0, whereas the latter limit 
agrees with Ref. The formal expansion in the scatter- 
ing matrix S, however, is only allowed if the probability 
of coherent (back)scattering from either the quantum dot 
or the contact is small, and does not describe interaction 
corrections to the conductance of a fully coherent quan- 
tum dot (i.e., a dot with a unitary scattering matrix). 
With the same arguments, the result for the interaction 
correction to the dot's capacitance, which was calculated 
in Ref. 's' using a formal expansion in S', is not applicable 
to a quantum dot with a unitary scattering matrix. The 
correct theory, in the limit of large iV, will be given here. 
We will not compare non-perturbative results in the two 
formalisms, that were obtained using refermionization^ 
or instanton solutionsi2ii2^ 

With regard to the interaction corrections to conduc- 
tance, our result is particularly simple: for a fully coher- 
ent dot there is no correction to the ensemble averaged 
conductance, to leading and sub-leading order in 1/iV. In 
particular, this applies to weak localization, the quantum 
correction to the ensemble averaged conductance that is 
suppressed by the application of a magnetic field. The 
reason for the absence of an interaction correction to the 
ensemble averaged conductance is the following. In the 
Landauer description of transport, conductance is related 
to a sum of squares of moduli of elements Sij (e) of the 
scattering matrix S{e), appropriately integrated over en- 
ergy if the temperature is nonzero (see, e.g., Ref. i4tf l. 
The ensemble average of the conductance depends on the 
ensemble averages (|S'(e)ijp) only. According to the ran- 
dom matrix theory of quantum transport, such an aver- 
age is determined by symmetries only (presence or ab- 
sence of time-reversal and spin-rotation symmetry). At 
leading and sub-leading order in 1/A'', we find that the 
interaction corrections can be expressed through a renor- 
mahzed scattering matrix: S{e) S'{e) = S{s) + 6S{e) 
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[see Eq. (|3.6|l below]. Although this interaction correc- 
tion changes the scattering matrix, and hence conduc- 
tance, of a specific dot, it does not change the symmetry 
of the scattering matrix ensemble, leaving the ensemble 
averaged conductance (including the weak-localization 
correction) unchanged. We will revisit this argument in 
more detail in Sec. IVIll 

The outline of our paper is as follows. In Sec. |n] 
we present a detailed discussion of the relevant exist- 
ing works in the literature. This discussion will serve to 
introduce the necessary concepts and notations, and to 
explain the disagreement between the two existing the- 
ories of 'weak Coulomb blockade' in coherent quantum 
dots. In Sec. IIIII we proceed with a precise definition of 
the problem and an exposition of the formalism. In the 
diagrammatic language that we will introduce later, the 
leading interaction correction to the current arises from a 
Fock-type diagram. The idea that interaction corrections 
to conductance in the presence of impurities should be 
thought of as due to scattering from an Hartree-Fock po- 
tential was introduced in the work of Matveev, Yue, and 
Glazmani^ii^ In order to make contact to that work, we 
have added to Sec. IIIII a calculation similar in spirit to 
that of Refs. l41i.42i that already contains the basic struc- 
ture of the interaction correction to the current. The 
full calculation of the current through the quantum dot, 
for the general non-equilibrium and time-dependent sit- 
uation is then described in Sec. IIVI Sections and IVll 
contain a detailed analysis of these results for nonlin- 
ear steady-state transport and adiabatic time-dependent 
transport, respectively. We conclude in Sec. I VIII Finally, 
appendices ^ El ^^nd El contain materials on the corre- 
lation functions in the 'one-dimensional' formalism, the 
ensemble average, and the case of a quantum dot with 
partially coherent scattering, respectively. A summary 
of our results, focusing on the relation between our work 
and previous approaches to 'weak Coulomb blockade' ap- 
peared in Ref. |4^ 



II. PREVIOUS WORK 

In this section, we present a more detailed overview of 
the results published in the literature. This will also serve 
to introduce some of the needed concepts and notations. 

A schematic drawing of the system under consideration 
is shown in Fig.^ It consists of a quantum dot coupled to 
two electron reservoirs by means of point contacts. The 
point contacts have A^i and N2 channels each (counting 
spin degeneracy). The two contacts are characterized 
by energy-independent reflection matrices rd and rc2, of 
dimension A^i and N2, respectively. In this notation, the 
conductance of each contact is 



point 
contact 1 




point 
contact 2 



FIG. 1: Schematic drawing of the system under consider- 
ation: a quantum dot (center) coupled to source and drain 
reservoirs (left and right) via point contacts. The two point 
contacts have A'^i and propagating channels at the Fermi 
level. Direct reflection at the point contact is described by 
energy- independent reflection matrices rd and rc2, so that 
their dimensionless conductances are Gi = (e^/27rfi)(A*'i — 

tr''ci''Ii) G2 — {e^/h){N2 — trrj.2r^2)i respectively. In- 
teractions in the dot are described by the charging energy 

fie 



with 51 + 32 ^ 1 is referred to as 'closed'; an open dot 
has gi+g2 ^ 1- For future reference, we define the total 
number of channels in the point contacts N = Ni + N2, 
the total dimensionless conductance of the point contacts 
g = gi + 92, and the N x N combined reflection matrix 
of both contacts 



rci 
r,2 



(2.2) 



The discussion of this section will be limited to the 
linear dc conductance G of the quantum dot. Our dis- 
cussion of the literature will be limited to the cases of a 
'coherent' quantum dot, for which transport properties 
are expressed in terms of its iV-dimensional unitary scat- 
tering matrix S{e), and an 'incoherent' quantum dot, for 
which the dot itself is treated as a reservoir in equilib- 
rium, with a self-consistent chemical potential. ^'^ In the 
'incoherent' limit, coherent propagation across the dot, 
or backscattering from within, is not taken into account. 

Neglecting interactions for electrons on the dot, for the 
'incoherent' case the conductance is nothing but the clas- 
sical series conductance Gd = {e^ /h)gig2/ g of the two 
point contacts. For the 'coherent' case, the conductance 
is given by the Landauer formula, which we write as 



Gn — 



N1N2 

2^ N 

=2 



(2.3) 



2TTh 



de 



dm 

de 



trAS'(e)AS'1'(£ 



Here /(e) = 1/(1-1- exp(e/r)) is the Fermi distribution 
function and the matrix A is defined as 



i = 1,2. 



(2.1) 



znn 

A contact with Ni = gi, i = 1,2 is called ballistic or 
'ideal'; otherwise the contact is called 'nonideal'. A dot 



A,,, = ^ X 



N2 forj = l,...,7Vi, 
-Ni for j = 7Vi + l,...,iV, 



(2.4) 



The introduction of the matrix A is a convenient way to 
separate the classical conductance for ballistic contacts 
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[the first term in Eq. (|2.3|l ] from the quantum corrections 
and the effect of nonideal contacts [second term in Eq. 

The Landauer approach is an intrinsically single- 
particle picture, and as such cannot capture any effects 
of electron-electron interaction by itself. In a quantum 
dot, the most important part of the interaction al?!?^ is 
described by the Hamiltonian 



(2.5) 



where the charging energy Ec — e^/2C is related to the 
geometric capacitance C of the quantum dot, N^ot is the 
number of electrons on the quantum dot, and J\f = CVg/e 
is a dimensionless gate voltage used to set the equilibrium 
value of A'^tjot- 

The charging interaction (|2.5(l is responsible for the 
phenomenon of Coulomb blockade in closed dots at low 
temperatures: suppression of the conductance, except if 
the gate voltage Af is tuned to a point of charge de- 
generacy. As a dot is opened, a situation best realized 
in lateral semiconductor dots where point contacts may 
be controlled electrostatically, strong quantum charge 
fluctuations lead to the progressive diminishment of the 
Coulomb blockade. A theory of this 'weak Coulomb 
blockade' was developed by Flensbergji MatveevjS and 
Furusaki and Matveev)^ for the 'incoherent' case. For a 
quantum dot with two nearly ideal single channel point 
contacts (|rci|, |rc2| <C 1) and at temperatures T <C -Ec, 
the interactions were found to change the conductance of 
the contacts according to cjj gj + Sgj, j — 1,2, where 



,£(3/4) 

'r(i/4) 



J -1,2, (2.6) 



up to corrections of order Ird"*. Here C « 0.577 is the 
Euler constant. As in the non-interacting case, the con- 
ductance G of the quantum dot is simply the series con- 
ductance of the two point contacts, 



G 



e- 
2^ 



r(3/4) e<^E, 



r(l/4)V ttT 



-tr r^rl 



(2.7) 



The first term is the classical conductance with A^i = 
N2 = 2, accounting for spin degeneracy. The second term 
shows the renormalization of the backscattering at the 
contacts, see Eq. (|2.6|) . At low temperatures the system 
enters into a non-perturbative regime. We refer to Refs. 
I1I I3I4I I7I for the expressions for the interaction corrections 
at different numbers of channels in the point contacts. 

A different approach to the same problem was taken by 
Golubev and ZaikinS^, and by Levy-Yeyati et al.^^ who 



described the interaction Hamiltonian (|2.5() through an 
effective theory for the electromagnetic environment of 
the junction. These authors also calculated the interac- 
tion correction to the conductance of the point contacts 
for an arbitrary transparency of the contacts using 
as a small parameter. They found 



^9] 



3T ' 

g 27r2T ' 



T > gE, 



(2.8) 



where si and S2 are the so-called Fano factors of the 
contacts, 



9j 



J = 1,2. 



(2.9) 



The appearance of the Fano factors is quite common to 
corrections to the conductance that are perturbative in 
the interaction, see Ref. liHli^ In fact, Si is the only 
quadratic function of the product r^^rl^ that vanishes 
both in the limits Vd — > (ideal contacts) and rd — > 1 
(completely closed contacts) . Equation H2.8|l agrees with 
the large-iV generalization of Eq. H2.7|l for nearly ideal 



contacts 
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The results of Refs. Illl2l l3l4l being non-perturbative in 
Ec, are of conceptual importance, but describe the inco- 
herent case only. They neglect entirely the possibility of 
coherent backscattering from within the dot, or coher- 
ent transmission from one contact to another. Thus Eq. 
(|2.7|) misses quantum interference effects, such as weak 
localization (the small negative magnetic-field dependent 
quantum interference correction to the ensemble aver- 
aged conductance) and conductance fluctuations, which 
one expects to be described by the Landauer formula 
l|2.3fl if the conductance of the contacts is large enough. 
Interaction corrections to quantum interference correc- 
tions cannot be captured by a simple renormalization of 
the point contact conductances; one needs to consider 
the total conductance G of the quantum dot. Only at 
temperatures for which relaxation processes in the dot 
are dominant one expects the separate renormalizations 
of the point contact conductances of Eqs. 1)2. 6|l or l|2.8|) 
to be a sufficient description. 

The first results for the 'coherent' case, concerning in- 
teraction corrections to transport in open dots described 
by a unitary scattering matrix, were obtained by Brouwer 
and Aleineri (see also the review Ref. |^. They found 
G = Go + Gba with 
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Gba = 



{2a + T2+ ti)t:^T^ f sinh[(Ti + T2 + a + to)TrT] sinh[(fT - ta)nT] 



^ L '^'^sinh[(CT + ri)7rT]sinh[((7 + T2)7rr] 



sinh[(ri + to)'7rr] sinh[(r2 + tQ)7rT] 



l/N 



.(2.10) 



Here to ~ ir/EcNe'^ is a charge relaxation time and S{t) 
is the Fourier transform of the scattering matrix S{s), 
which is defined as 

S{r) = I ^S{e)e-^^\ (2.11) 

S\r) = J ^S\e)e-^^\ (2.12) 

For a dot with two single-channel spin-polarized point 
contacts, A^i = A^2 = 1, there is an additional interaction 
correction to the conductance that depends explicitly on 
the gate voltage N . Equation (|2.10|) is obtained as a 
formal expansion in the scattering matrix S. Such an 
expansion is controlled if S is subunitary, as is the case, 
e.g., in an effective description of inelastic processes in 
the dot 4^ Noting that only times up to the 'thermal time 
horizon' h/T contribute to the interaction correction of 
Eq. (|2.1U|) , Brouwer and Aleiner argued that their formal 
expansion is also justified for unitary S if the temperature 
is sufficiently high that the contribution from scattering 
times < h/T is small. For a chaotic dot with mean dwell 
time Td and ideal point contacts, this translates to the 
condition T ^ h/Td. 

The result of Ref . captures the Landauer part of the 
conductance exactly, whereas the interaction correction 
(|2.1l)|l reproduces the singular behavior of Eq. H2.7|l for 
the scattering matrix S{t) = rc(5(r) -f . . ., provided that 
the scattering from inside the dot, represented by 
is from far beyond the 'thermal time horizon'. Using 
known statistical properties of the scattering matrix for 
chaotic quantum dots, Brouwer and Aleiner calculated 
the ensemble averages of the conductance and the con- 
ductance fluctuations for T 3> h/T^, and found that both 
the weak localization correction to the conductance and 
the conductance fluctuations are slightly enhanced by in- 
teractions. 

In this paper we question the arguments used to jus- 
tify the formal expansion in the scattering matrix used 
to obtain Eq. (|2.1U|) . As we will argue below, the interac- 
tion correction of Eq. H2.1U|I is not the correct interaction 
correction to the conductance of a coherent quantum dot 
at any temperature, despite the fact that it is small if 
T ^ h/r^. Our main arguments will be given in the 
next sections, where we present two calculations of the 
interaction correction to the conductance: The first is 
a first-order-in- i?c calculation of the conductance, which 
differs from Eq. (|2.1U|) if that result is expanded in E^. 



The second is a calculation to all orders in Ec, but to 
leading order in 1/7V. This calculation confirms that Eq. 
(|2.1U|) is found if the amplitude for coherent scattering 
from the dot is small — i.e., transport through the dot is 
mainly incoherent (but not fully incoherent) — , but we 
find a different interaction correction if transport is fully 
coherent. 

The main shortcoming of Eq. (|2.1U|) can already be 
seen noting that, although for T ^ h/r^ most scattering 
processes have delay time ^ h/T and thus do not con- 
tribute to the interaction correction, the remaining in- 
teraction correction is dominated by scattering processes 
with delay time < h/T. Hence, it is necessary that the 
theory describes this time range accurately; its contri- 
bution being small is not sufficient to justify Eq. H2.10(l . 
That Eq. H2.10|l does not describe delay times < h/T 
accurately can be seen by considering the example of a 
'dot' that consists of a single ballistic channel in which 
there is an interaction of the form H2.5|l . The channel 
is of length L, so the 2x2 scattering matrix describing 
ballistic propagation from one end to the other is 



S{r) 



1 

1 



(2.13) 



where Td = L/vp, vp being the Fermi velocity inside the 
channel, is the time to propagate through the channel. 
We assume Td » to. Then, substituting Eq. (|2.13|) into 
(|2.1UI) gives a conductance that is larger than the con- 
ductance quantum e^/27rfi,, G — e'^/2TTh + Gba, with 



Gba 
Gba 



TrTTd 



2TTh h 

irh 

2^ STrd 



g-2.Twr.^ T^h/Td, 



T<t:h/Td. (2.14) 



On the other hand it is well known that in the absence 
of backscattering the conductance of such a system can- 
not deviate from the quantized value /h, see, e.g., Refs. 
45-46 41 so that the interaction correction to the conduc- 
tance must be zero for all temperatures. Hence, although 
the error made by using Eq. H2.1()|l is exponentially small 
for Td ^ h/T , it is substantial for Td < h/T . 

Very recently, Golubev and Zaikin extended their envi- 
ronmental formalism to the case of a fully coherent dot.^'* 
Again using 1/A^ as an expansion parameter, they found 
G Go -I- Ggz, with 
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GZ 



Im / dedio 



df{e 



K{uj)tr [AS'(e)AS'^(£ + w) - AS{e)S\e + u;)S'(£)AS'^(e)] (1 - 2/(e + u;)), 

(2.15) 

I 



where k(w) is an effective interaction kernel that de- 
scribes the fluctuations of the dot potential. Since these 
fluctuations are small as 1/A^ if N is large, the correc- 
tion H2.15|l is a small correction to the Landauer result 
Go . Golubev and Zaikin performed the ensemble average 
over S and found that the leading interaction correction 
for the coherent case is the same as in the incoherent 
limit (i.e., the classical combination of the corrections 
Eq. \2.ii\\ ) as long as T ^ h/T^, and saturates when 
T < ft/Tfj. The agreement with the incoherent limit is 
no surprise, since the leading-in-A'^ contributions to the 
conductance are generally insensitive to the presence of 
quantum coherence. 

The authors of Ref. 24 did not analyze the quantum 
interference corrections to the conductance and thus at- 
tempted no comparison with the results of Refs. 00. 
Our full calculation, which is described in the next sec- 
tions and which includes quantum interference correc- 
tions, agrees with Eq. (|2.15|) and extends it to the case 
of general time dependent transport — the focus of this 
paper. 

Technically, our formalism is very close to that of Refs. 
Within the same formalism, we can obtain Eq. (|2.1U|) 
using a formal expansion in the scattering matrix S and 
Eq. H2.15|l as a formal expansion in 1/A^. As we dis- 
cussed above, expansion in powers of S requires a relax- 
ation mechanism, and Eq. (|2.1U|) can be justified only if 
relaxation is strong, so that, effectively, \S{e)\ <C 1. 



III. FORMALISM 

A. Definition of the problem 

The system under consideration — a quantum dot cou- 
pled to source and drain reservoirs via point contacts — 
has been introduced at the beginning of Sec. ^ There 
are two point contacts, with A'^i and N2 channels each (in- 
cluding spin degeneracy). Electrons on the dot interact 
via the sim ple interaction Hamiltonian H2.5|l . We refer 
to Refs. 35 for a microscopic justification of Eq. (|2.5|) 
for quantum dots with a large (internal) dimensionless 
conductance. In using Eq. (|2.5|) we ignore those parts of 
the interaction describing superconducting pairing and 
exchange]^ 

Following Flensbergi and Matvee\i^ we describe the 
electron dynamics in the point contacts by a one- 
dimensional Hamiltonian. We locate the lead-dot inter- 
face at a; = 0, taking dot and lead at a; > and x < 
respectively, see Fig. |21 Since there is no backscattering 







lead 1 
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lead 2 
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quantum 
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FIG. 2: Representation of the source and drain reservoirs 
as one-dimensional ideal leads. The lead-dot interface is at 
2; = 0. The one-dimensional description is valid from 

a; = — cx) up to slightly beyond a; = 0. 



of electrons that have left the dot through the point con- 
tacts and into the reservoirs back into the dot, we may 
represent the reservoirs by extending the one-dimensional 
description to all x < 0, 



H = ivp 



N 



dx 



d 



V']l(^)^V',l(^) (3.1) 



(x) 



Here the index j labels the channels in the two point con- 
tacts, including spin, TV is the total number of channels 
in the two point contacts, ipjL and ipjfi are annihilation 
operators for left-moving and right-moving electrons in 
channel j, respectively. Because of scattering by the dot, 
these are not independent fields. In the one-dimensional 
Hamiltonian (|3.1|l . the kinetic energy is linearized around 
the Fermi energy, which is appropriate if the Fermi en- 
ergy is not too close to a threshold at which a new channel 
in the point contact is opened. The coordinate x in the 
one-dimensional Hamiltonian H3.1|l arises as the Fourier 
transform of the one-dimensional Hamiltonian in momen- 
tum representation, see Ref. 1^ Except for the immediate 
vicinity of the point contact, it cannot be identified with 
a distance to the lead-dot interface. 

It is convenient to follow Aleiner and Glazman,- and 
write the number of electrons on the dot as 



dot 



E 



dx 



V']L(2;)V'jL(a;) 



(3.2) 



Here iVref is the total number of electrons in the system 
(reservoirs and dot), which is time independent. Now the 
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interaction does not contain the dot degrees of freedom. 
Thus the relation between the operators for right mov- 
ing electrons entering the quantum dot and left moving 
electrons exiting the quantum dot is given by the non- 
interacting formula. In the general time-dependent case 
one then has^ 



N 

k=l •' 

N 



fe=l • 



(3.3) 



where S{t, t') is the scattering matrix of the dot. Notice 
that the scattering matrix has two time arguments. For 
time-independent transport, S{t, t') is a function of the 
time difference t — t' only. 

The relation 1)3. 3|l . together with the one-dimensional 
Hamiltonian (|3.1|) . Eq. (|3.2|) for N^ot, and the bound- 
ary conditions set by the reservoirs completely define the 
problem. Notice that we never explicitly introduce the 
dot degrees of freedom into the problem. If the scatter- 
ing from the dot is not fully coherent, Eq. (|3.3(l has to 
be replaced by a different boundary condition. We re- 
turn to the case of incoherent scattering in section UlI Dl 
A model for partially coherent scattering is discussed in 
appendix lO 



B. First order in Ec 

It was shown by Matveev, Yue, and Glazman^ii^ that, 
at first order in the interaction, interaction corrections 
to conductance in the presence of impurities may be un- 
derstood in terms of the elastic scattering of electrons 
from the Hartree-Fock potential created by the impu- 
rities. This gave a simple and intuitive picture of the 
renormalization of backscattering in the interacting one- 
dimensional electron gas, without recourse to sophisti- 
cated formalism. We now outline a similar calculation 
for transport through a quantum dot. 

For linear transport, it is sufficient to consider the 
scattering off the Hartree-Fock potential in equilibrium. 
With the replacement (|3.2I) the Hartree-Fock potential 
acts in the leads, not in the dot. In equilibrium there 
is no Hartree potential as (N^ot) — A/". For the Fock 
potential we find 



N 



Hp = 



/ dxdy '0|^(x)V'y(a;,?;)'0^.^(y) +h.( 



(3.4) 

where "h.c." denotes the hermitian conjugate, plus for- 
ward scattering terms that do not affect the linear con- 
ductance at lowest order in Ec- The Fock potentials are 
written in terms of the density matrices, which may be 



evaluated using the boundary conditions (|3.3|) 

V,jix,y) = -2ES]r{v)Al{^)) 

Jo 

iT 

X 



2sinh [ttT {t - vf[x + y] + iX)] 



,(3.5) 



where A is a positive infinitesimal. Thus the effect of 
the Coulomb interaction is to establish a non-local Fock 
potential in the leads within a region of order the thermal 
length vp/T from the contacts. It is now straightforward 
to find the change in the scattering matrix due to this 
potential 

SS{e) = ^ J ^K{io)[{2f{e-uj)-l)S{e~uj) 

+ (2/(e -I- c^) - l)S{e)S\e + co)Sie)] , (3.6) 



where 



Er 



(3.7) 



0"*" being a positive infinitesimal. The first term in H3.6() 
represents the effect of direct i? — > L (in to out) scatter- 
ing by the Fock potential H3.4|l . The second corresponds 
to the dot scattering R L, followed by the Fock poten- 
tial scattering L ^ R, and ending with the dot scattering 
R ^ L once more. Substitution of 1)3.6(1 into the Lan- 
dauer formula H2.3|l immediate yields the result (|2.15|) 
for the interaction correction to conductance, with k,{uj) 
given by Eq. 1)3. 7|l above. 

We may repeat the deriv ation for the incoherent case 
that is considered in Refs.'l"?3"i'25"23. In this case, the 
inclusion of the Fock potential leads to a change of the 
reflection matrix for direct reflection at the contacts. 
Since electrons leaving the dot are incoherent with those 
entering it, they do not contribute to the Fock potential 
()3.5)) . The renormalized reflection matrix is thus 



<5rc(e) - ^ J ^^(0;) [{2/(6 ~ cj) ~ l)r, 
+ i2f{e + Lu)-l)r,M , 



(3.8) 



Taking into account the energy dependence of the renor- 
malized reflection matrix r^, we find that the interaction 
correction to the dimensionless conductance of contact j, 
j — 1,2, becomes 



2sjgjlm I de I ^ 



(3.9) 



dm 

de 



il-2f{e + Lu))4io), 



where sj was defined in Eq. 1)2.9)) . Although this looks 
very similar to 1)2.15)1 with the scattering matrix chosen 
as S{e) = Tc, the effect of using a non-unitary scattering 
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matrix is significant. Notice that for a unitary scatter- 
ing matrix, the trace in Eq. (|2.15|) vanishes as w^, so 
that the regulator in (j^.Tfl plays no role. In H3.9() on the 
other hand, the imaginary part of the integrand comes 
from the region near lu ~ and yields a finite, regulator 
independent contribution to the conductance 



^9]^-Sjg~, i = l,2, (3.10) 

which coincides with the high temperature limit of Eq. 
(|2.8|) .^» Recovering this same result, together with the 
low-temperature limit of (|2.8|) . as the zero level spacing 
limit of scattering from a coherent dot is subtle, see Sec. 
Eland Ref. Hi 

With this physical picture of the interaction correc- 
tions, we now proceed with the formal development of 
the theory. 



C. Effective action 

We are interested in the current through the point con- 
tacts, at the dot-lead interface x = {). The current in 
channel j is 



-evp 



^]fl(0)V^,fl,(0)-^JjO)V',x(0) • (3.11) 



0L\ 



Below, the expectation value of the current is calcu- 
lated using standard methods of nonequilibrium many- 
body theory,''^ The current I{t) is related to the state of 
the system at a reference time tref at which the Hamilto- 
nian is not time-dependent. The expectation value I{t) is 
then found from a thermodynamical average at the time 
t-cei- An important issue here is that the total number of 
electrons, N^ci, be kept constant — especially when ob- 
servables at two different values of the gate voltage N are 
compared,^ In Ref. this was implemented by taking this 
thermodynamical average in a canonical ensemble. In our 
case, however, we can take the average in a grand canon- 
ical ensemble, since, in the non-equilibrium formalism, it 
is sufficient that the dynamics conserves the number of 
electrons. 

The one-dimensional dynamics of the Hamiltonian 
(|3.1|) applies to the region x < only; the dynamics in 
the quantum dot is described by the boundary condition 
(|3.3() . This boundary condition is difficult to implement 
in the calculations, and it is advantageous to reformulate 
the problem in terms of a one dimensional Hamiltonian 
of the form H3.1|) . where the x integration extends over 
the entire real axis. In this formulation of the problem, 
the boundary condition H3.3|l . or its equivalent for a dot 
in which scattering is not fully coherent, is replaced by an 
effective action iS that involves the fermion operators at 
X = rj, where 77 is a positive infinitesimal. Such a formu- 
lation of the problem has first been given by Aleiner and 
Glazman,^ see also Ref. 8, Here we will use an equivalent 
but simpler effective action formulation. 



In order to derive the effective action, a fictitious half- 
infinite A^-channel lead is side-coupled to the point con- 
tacts, at position x = t], where 77 is a positive infinitesi- 
mal. The fictitious lead couples only to electrons moving 
out of the dot, i.e., to left movers. The total Hamiltonian 
H is written as 



(3.12) 



where Hq is the Hamiltonian of the full interacting sys- 
tem (including the quantum dot, the reservoirs, and the 
fictitious lead) , with ideal coupling between the fictitious 
lead and the point contact, see Fig. |3| (solid lines). Hence, 
in the system described by Hamiltonian Hq , all electrons 
entering the point contact from the dot will exit through 
the fictitious lead, whereas all electrons entering from the 
fictitious lead will exit towards the source or drain reser- 
voirs. The Hamiltonian Hi represents an impurity in the 
contact to the fictitious lead which ensures that electrons 
coming from the dot in channel j have amplitude —irj 
to be transmitted coherently towards the source or drain 
reservoirs (cf. dashed lines in Fig. |3J) and that electrons 
coming in from the fictitious lead have amplitude ~irj 
to be reflected into the fictitious lead, 

N 

Hi ^Y.'^VFVj [i^l^m^M + ^1^(277)7^^.^(0)] , 

(3.13) 

where 



27;,- 



j^l,...,N. 



(3.14) 



The parameters Vj (or Vj ) describe how well the fictitious 
lead is coupled to the point contact. In the limit — )■ 1, 
j — 1, . . . ,N, the fictitious lead is fully decoupled. We 
will take the limit Vj 1 at the end of the calculation. 
In writing Eq. (|3.13() , the same regularization of fermion 
operators has been employed as in Ref. 's'. 

The operator identity (j3.3|l applies to left-moving 
fermions at the moment when they exit the dot. With 
the inclusion of the fictitious lead, this means that one 
should take ipL sd, x — 277, not at a: = 0. Hence Eq. 1)3. 3|l 
gives a relation between 7/1^(277) and 7/;fl(0), 

^ r 

k=l •' 

(3.15) 

The advantage of separating the Hamiltonian into the 
contributions H^ and Hi is that the problem described 
by the Hamiltonian Hq alone is exactly solvable, see Sec. 
nil El The effect of Hi is then treated in perturbation 
theory. Using Eq. (|3.15|l . the expectation value of the 
current at time t can be expressed in terms of an effective 
action S, 



(3.16) 



/,(i) = (Tce-''^/,(i) 
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reservoirs 



x=0 T| 2t| 



ref 



reservoirs R 




Y fictitious lead 



R 



L 




V fictitious lead 



FIG. 3: Inclusion of a fictitious lead at the point contact. Top 
panel: With the Hamiltonian Hq, the fictitious lead couples 
ideally to outgoing electrons (solid lines), so that all electrons 
that exit the dot enter into the fictitious lead. The Hamil- 
tonian Hi describes scattering of electrons exiting the dot 
into the physical lead and backscattering of electrons coming 
from the fictitious lead (dotted line). The amplitude for this 
scattering process in channel j is irj, j = 1, . . . , A^. Bottom 
panel: If rj = l, the fictitious lead is decoupled again and the 
original problem is restored. 



FIG. 4: Integration contour for Eq. 



to infinity and set 

N 

Ho = ivF 



dx 



(3.18) 



where, as before, iVdot is expressed in terms of the fermion 
fields for X < only, see Eq. H3.2() . In the system de- 
scribed by the Hamiltonian Hq, the fermion fields ipu 
have chemical potentials fiR corresponding to the chemi- 
cal potentials in the reservoirs, whereas the fermion fields 
ijjL have the chemical potential /ii of the fictitious lead. 
For the physically relevant case Vj ^ 1, j = 1, . . . , N, the 
fictitious lead is fully decoupled from the point contact. 
In that limit, all physical observables become indepen- 
dent of the choice of fiL ■ 

The effective action j[3T7|l is different from the effective 
action used in Refs. lyQS the kernel in the latter action 
is a matrix L, which is a function of the scattering matrix 
S of the quantum dot. 



where the current operator Ij is given by Eq. (|3.11|) 
above, "c" denotes the Keldysh contour, see Fig. 0] Tc 
denotes contour ordering along the Keldysh contour, the 
time-dependence of the operators is that of the interac- 
tion picture with respect to Hq, and the average (. . .)o 
denotes an average with respect to Hq at reference time 
trei- The effective action S is 



S 



(3.17) 



N N 



2vF / dti / dr^^v 

•''^ '' j=l k=l 

X V-Il (0, ii ) 1 , ti - T)?Afcfl (0, i 1 - t) 
+ 4i?(0, h " r){S^)k, ih - T, ti)^,x(0, ti) 



where we dropped reference to the infinitesimal 77. It 
is important to note that in this effective action, the 
contour-ordering occurs according to the 'contour time' 
ii, not according to the 'scattering delay time' r. 

Both the current operator (|3.11l) and the effective ac- 
tion H3.17|l do not contain any reference to the electrons 
after they have passed through the point contact, so that 
one may send the upper integration boundary in Eq. (|3.1I) 



D. Relaxation inside quantum dot 

If the fictitious lead used for the above derivation of 
the effective action is fully coupled to the point contact 
{i.e.,, if = 0, j = 1, . . . , N), every electron exiting the 
quantum dot will escape towards the fictitious lead, not 
towards the (physical) reservoirs. Setting the chemical 
potential /il of the fictitious lead such that no net current 
is drawn, one finds that each electron that escapes into 
the fictitious lead is replaced by another one — without 
phase relationship. Hence, coupling to the fictitious 
lead allows for complete phase and energy relaxation of 
electrons exiting the dot, just before they pass through 
the point contact to the reservoirs. 

The simplest model for incomplete relaxation in the 
quantum dot is to use the coupling to the fictitious lead 
with < Vj < 1. In this case, the electrons are al- 
lowed to escape into the fictitious lead and relax, but 
with finite probability only. Hence imperfect coupling to 
the fictitious lead describes a quantum dot with a finite 
probability that the electron exits the dot without relax- 
ation. Of course, this model is over-simplified, because 
all relaxation is localized at the point contact. A more 
realistic, but still phenomenological model for relaxation 
in the quantum dot is discussed in appendix [HI 

A connection can be made between the effective action 
(|3.17(l for general Vj and the 'fully incoherent' limit of 
the problem, with nonideal contactSfiaMiMM? Hereto, 
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the fictitious lead is interpreted as the "quantum dot', 
the scattering matrix Sij(t, t') is replaced by the identity 
matrix, Sij{t,t') — SijS{t — i'), and one identifies the 
reflection matrix r with the contact's reflection matrix 
Tc. In this case, energy and phase relaxation inside the 
dot is still complete, but a nonzero value of the reflection 
probabilities r| allows one to describe a nonideal coupling 
between dot and reservoirs. 



For right movers, however, the presence of the charging 
interaction in the lead means that the ^electrochemical 
potential fic is the relevant quantity 



(/5,fl(0,t)> 
Aicj(i) 



1 



-Mci(i), 



2TrvF 



(3.20) 
(3.21) 



E. Correlation functions 

In order to formulate a perturbation theory in the effec- 
tive action (|3.17|l . we need the contour-ordered fermion 
correlation functions of the Hamiltonian Hq- The equilib- 
rium correlators were calculated using the bosonization 
method in Ref. 1^. The new elements in the nonequi- 
librium case are few. We discuss them now in physical 
terms, with the details being relegated to appendix IXI 

We specify the non-equilibrium conditions applied to 
the system through the chemical potentials fijLit+x/vp) 
and fj,jR(t — x/vp) of the left and right moving electrons. 
Although the chemical potentials of the flctitious leads 
should not appear in physical quantities in the Vj — > 1 
limit, we retain them in order to illustrate how this hap- 
pens. 

In the absence of scattering, i.e. setting Hi = 0, the 
Hamiltonian is quadratic and a mean-field treatment is 
thus exact. The average density of left moving electrons 
at the contact is simply related to their chemical poten- 
tial 



1 



2TrvF 



(3.19) 



where Vd (t) was defined in (|1.1() . The average current at 
the contact is 



I,{t) - ^evF{{p,BiO,t))-{pjL{0,t))) 



(3.22) 



e 

2^ 



fijR - PjL - ((iVdot(i)> - AA) 



Thus with an excess charge in the dot (/Vdot) A/", there 
is a contribution from the interaction. We solve (|3.20|l 
by writing the excess charge as an integral of the current 



2ttvf 



(3.23) 



TT ^ 

fc=l 



dt'{pkR{t')) - {Pkdt')), 



where we have introduced the convention that fields with- 
out a position label are taken to be evaluated at the ori- 
gin. The solution of (|3.23() is 



2ttvf 



PjL 



(i)+V' / drijkir) {pkiiit - t) - pkL{t - t)) 



(3.24) 



where we defined the kernel 



imkir) = SmkS{T - 0+) 



0{r), (3.25) 



0+ being a positive infinitesimal. The kernel ijkir), 
which appears frequently in the following development, 
relates the current at the contacts to the chemical po- 
tentials in the leads, accounting properly for the effect of 
the charging interaction. Note the appearance of the RC 
time tt/EcN characteristic of an ideal A^-channel contact. 

Combining the above results, we find that the current 



in channel j is given by 



N 



X [fikRit - t) - HkL{t - t)] . (3.26) 
Note that in the dc bias situation Eq. ()3.2tif) simplifies to 



e 

'2^ 



N 



fc=i 



1 

N 



{P'kR — Pkh) , 



SO that for the two-terminal set up 



e 

2^ 



(3.27) 



(3.28) 
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where A^i_2 = ^J■Rl,2 — Mil, 2 are the chemical potential where A is a positive infinitesimal and 
differences in the two contacts, and is the matrix 
defined in Eq. (1231). 

Correlators of many fermion operators will be ex- 
pressed in terms of the single- fermion Green functions, 

-*(T^c^™L(s)^iLW) = GmnL{s,t) (3.29) 
-*(7^c^™flW^L(5)> - GmnR{t,s). (3.30) 

These read 

GmnL{t,s) : 



2vf sinh[7rr(t — s — i\sign^{t — s))] ' 

A 7^p-i(0mH(t)-0mn(s)) 

Q ^-j "rnn-'- 

' 2vp sinh[TTT{t — s ~ iXsigii^{t — s))]' 

(3.31) 



Pm.L{t) = 
t'm.Ri.t) = 



dT^rnL(t - t). 



(3.32) 



/ dr fimL{t-T)+y2 dT'i,nk{T'){HkR{t - T - t') - HkL{t- T - t')) 

Jo I ^^^Jo 



are the integrals of the electrochemical potentials of right and left movers. It is straightforward to see that (|3.3HI 
is consistent with (|3.19(l and l|3.24(l . The single- fermion Green function can be written as a sum of a Keldysh part, 
which does not depend on the contour positions, and a delta function. 



Ik * 

GmnL{t,s) — -^G^^];^{t,s) — -^—6mnSc{t,s), 

1 K * 
GmnR{t,s) — —G^^j^{t,s)~- SmnSc{t,s). 

2 ZvF 

Here we abbreviated 

Sc{s,t) = sign^{s - t)6{s - t). 
Now the contour-ordered average of a product of fermion operators is calculated as 



(3.33) 



(3.34) 



P,Q i=\ j = l ij 



l/N 



(3.35) 



where P and Q are permutations of the numbers z = 1, . . . , n and j = 1, . . . , m, respectively, and the function f{t — t') 



IS 



\nf{t-t') = In 



XEcNe^ EcN 



-BcivcA, sinh[7rT(t - t' + ( - zO+sign,(t - t'))] 
^ sinh(^rC) 



(3.36) 



In the last equation, 0+ is a positive infinitesimal and C « 0.577 is the Euler constant. We separate ln/(t — t') into 
a contour- independent (Keldysh) part and a contour-dependent part. 



\nf{t - t') = In |/(t -<')!+ ^NKo{t - t')sign^{t - t'), 



(3.37) 
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where 

Ko{t -t') = ^[l- e-^^^^(*'-*)/-) e{t' - t). (3.38) 

Here 9{x) = 1 if x > and otherwise. 

Note that the result (|3.35|l is compatible with the symmetry of the Hamiltonian Hq. At first glance this is U{N) x 
U{N), the two factors corresponding to the right and left movers. Certainly this is the case for the non-interacting 
model with Ec 0, where (|3.35|l reduces to the usual Wick's theorem result. The interaction, however, leads to 
the chiral anomaly, as in the Schwinger model, that reduces the symmetry to SU{N) x SU{N) x C/(l). Physically, 
this is because the interaction can scatter right movers to left movers and vice versa. This allows a second type of 
correlation function to exist, characterized by the 5?7(A^)-invariant antisymmetric tensor £i-^i2...iK- We will refer to 
such correlators as anomalous. They describe the scattering of right movers to left movers (and vice versa) by the 
interactioni^ Below we list the first of these nontrivial correlators for = 1 and for N = 2. For N — 1 one has 

-^{TJ,mU-')) = -:|gt£0e-(^-W-^-(^')-2'^^). (3.39) 
For N = 2 the relevant anomalous correlators are 

(-*2) (Tc^iaXO, ^1)^-^(0, i2)4fl(0, i3)^2L(0, t4)] 

XT' C \ ^ 



i-i^) (TcViiJO, S,)lP\j,{0, 4)V'2l(0, 53)V'L(0, 4) 
Ere 



^ ^^iMs^)-4>Ms',)+4>L.is.)~Ms',)^2n^r (^^^^ s',) f (s,, s',)f{s,, 4)/(s3, s',)f^ . (3.40) 



Notice that the 'vacuum-angle' describing the U{1) sym- 
metry breaking is just 2TrM. 

For the perturbation theory calculation of the current, 
we need contour-ordered correlators of the form 

{Tj,{t)A{t[)...A{t'J), 



respectively. Factors Em for the creation operators for 
left-moving and right-moving fermions in channel k are 
-Fj^kL and -Fj^kR, respectively. 



where the symbols A{t') represent creation or annihila- 
tion operators for left or right moving fermions. There is 
a simple relation between such correlators and the corre- 
sponding correlator without current operator, 

{Tj,{t)A,{t[)...AM) 



n 



(Teii(t'i)...i„(C)) J2 ^™(^'C)- (3.41) 

m—l 



Here the Fm depend on whether the corresponding oper- 
ator Am is a creation or annihilation operator for left or 
right moving fermions. For an annihilation operator for 
left-moving and right-moving fermions in channel k, one 
has 

Fj.khitM) = -— / dTijfe(t-ti-T)coth(7rTT) 



-sign^(t - ti)ijk{t - h) 



(3.42) 



IV. CALCULATION OF THE CURRENT 



With the help of the correlators listed in the previous 
section, the perturbation expansion (|3.1()l) can be evalu- 
ated. In this section, we first describe a full evaluation of 
I{t) up to second order in the effective action (first order 
for A^ = 1). We then reorganize the perturbation theory, 
and calculate I{t) to all orders in 5, but in a perturba- 
tion expansion in 1/A^, which turns out to correspond to 
a loop expansion!^ 
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A. Expansion in effective action 



To first order in the effective action one has 



To zeroth order in the effective action, one finds that 
the current in channel j is given by 



(4.1) 



e 

2TTh 



N 

/ dh ^ ijk{t - ti) [nkL{ti) - ^J-kR{tl)] , 

k=l 



I J At) 



(4.2) 



Substituting Eq. H3.17(l for the effective action and using 
the fermion correlation functions of the previous section, 
in agreement with Eq. H3.26|l and with ijk{t~ ti) defined one finds Ij,i{t) = 0, except for = 1. In that case, one 
in Eq. (|3.25(l above. has 



hit) = AvFvlm J^dt, J dTe-^^'^S{tuti - T)[-FL{tM) + FuAM - r)] ^nCFfiO) ) e"'^^"^*^""^"^"^*^"- 

(4.3) 

In this equation, it is important to note that all contour-dependence should depend on ti only; the time difference r 
should not enter into the contour signs. We also note that the argument of /(r) is always positive, so that /(r) does 
not depend on the contour position, the only the contour-dependent part of the function F that contributes to the 
integration, and that the contribution from right moving fermions is zero by causality, so that 

hit) = -^M" / '^^'^^^'^'^ / (4-4) 

One verifies that for a time-independent situation — Sit,t — r) independent of t and chemical potentials fiL and fin 
independent of time — the current is zero. 

The second order calculation proceeds similarly. We first restrict our attention to the case N > 2. Using matrix 
notation, we then find 



Ij,2it) = --i2vF)^ j dtidsi J dridai 



/(n)/(ai) Y^'' 



Jiti- si+ai)fisi-ti+Ti), 
X trv'^[ijit, ii)sign^(i — ti) — ijit, Si)sign(,(t — Si)] 

X GLisi,ti)Siti,ti - Ti)GRiti - Ti,si - ai)S\si~ai,si). (4.5) 
Here ij is an N x N diagonal matrix, closely related to the kernel defined in Eq. (|3.25|) . 

i'^j)7iini'^) — ^mn'^jmi^)- (^■^) 

At this point, we write the Green functions Gn and Gl as the sum of the Keldysh component and a delta function, 
see Eq. (|3.33ll . Doing this, there will be four terms: one term with delta functions for both Gl and Gr, two terms 
with one delta function and one Keldysh Green function, and one term with two Keldysh Green functions. The first 
term is easily shown to vanish. In the second and third term, the interaction function / cancels, and one finds 



Ij,2,oit) = 2ivFe / dtitiv ijLit,ti) 



GlitiM) - j dTidaiSituh - Ti)Gg(ti - ri, ii - <Ji)S\ti - <JiM) 



(4.7) 



Here the last index "0" indicates that, apart from the difference between ijfe(T) and SjkSir), Eq. H4.7|l — together 
with the zeroth-order result of Eq. I|4.1|l — reproduces precisely what one would obtain from the scattering approach 
for non-interacting fermions. Finally, in the remaining term, the interaction contribution is retained. In fact, it is just 
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because of the contour-dependence of the function / that the one finds a nonzero result. We then find 



it) 



'2eivp J dti j dsi j dridai 



l/N 



/(ti-.si+ai)/(si-<i+Ti) 
X trv^ij{t,h) [G^{si,h)S{tuh - Ti)G^{h - n, si - (Ti)(5^)(si - cti, si) 
- Gf{h,si)S{si, si - ai)Gl{si -a^M' ri){S^){ti - t^M)] ■ 

For N = 2 there is an additional oscillating contribution to the current, 



sin[Ko(si - ti+ Ti)] 



(4.8) 



.it) 



ieviV2 



^2/(0)j 



Ime 



-2TTiJ\f 



dti / dt2dt[dt2sm{(t){t2,t[)) 



X det[zj(i,ii)e"^^(*i)S'(ii,t'i)e-*'^«(*'i\e*'^^(*^)5(t2,t2)e"*"^"'*^^] 

X |/(^l,^;)/(^l,^^)/(^2,t'l)/(^2,^i)^/^ 



where we defined 



det(i3(l),...,B(7V)) = E (-l)''S(l)fex,P(/cO---5Wfc«,P 



(few) ' 



P fci,...,fcjv 



(4.9) 



(4.10) 



where B{j) is an N x N matrix, j = 1, . . . , N. Equations 
(|4.4|) . (|4.8|) . and (|4.9|) represent the interaction correc- 
tions to the current. 

This result is found to agree precisely with the linear 
dc conductance (|2.1U|I calculated in Refs.0and0in the 
dc bias limit, to first order in the bias. To see this, we 
set HHi — fj,R2 = —eV, expand to linear order in the bias 
voltage, take S{t, t — r) to depend on t only, set fijL — 
for all channels j = 1, . . . , in the fictitious lead, use 
the approximations |/(t)| w AvrT/uf^ sinh[7rr|r -|- to|]i 
ko(t) = ttO{—t — to)/N valid for r ^ tg and r :§> to, 
and set Vj — 1, so that rj = 1 (see Eg. I3.14|l . which cor- 
responds to detaching the fictitious lead and thus restor- 
ing the original system. The oscillating current correc- 
tions H4.4|l and (|4.9|) reproduce the periodic-in-A/" interac- 
tion corrections to the dot's 'electrochemical capacitance' 
= dQ/dV found in Ref. i if we set fj,jL = 0, choose 
a time-dependent chemical potential fijuit) = eVt for 
all channels, and expand in V . [In order to recover the 
non-periodic interaction correction to the capacitance re- 
ported in Ref. 8 from Eqs. (|4.7|l and (|4.8|) one has to 
include time-dependencies in both and hr in order 
to ensure that no charge escapes through the fictitious 
lead.] 

However, truncating at second order in the effective 
action S is not a satisfactory description for a fully co- 
herent quantum dot as there is no a priori reason why the 
effective action is small in this case. This will be made 
explicit in the next subsection, where we report a calcu- 
lation to all orders in S (but for N large). However, that 
Eqs. H4.4I) . (|4.8I) . and (|4.9|) cannot be the correct inter- 
action correction for a fully coherent dot becomes clear 
when one realizes that the current in fact depends on the 
chemical potentials /ijx, j = 1,...,N, of the fictitious 
reservoir, despite the fact that one has taken the limit 



rj—^l,j — l,...,N. (Note that this dependence on the 
fictitious chemical potentials is no t m anifest in the lin- 
ear response calculation of Refs. I ^TlSl because fijL — 0, 
j = 1, . . . , A^, in equilibrium.) 

Of course, truncating the perturbation expansion at 
second order in S is justified if there is a reason why 
S is small. This is the case, e.g., if the rj are small, 
j = 1, . . . ,N. Then the fictitious lead is strongly coupled 
to the point contact, and serves as a source of relaxation. 
In that case the chemical potentials fijL must be chosen 
such that no current flows through the fictitious lead at 
any point in time. A more realistic (but still phenomeno- 
logical) model of relaxation is to couple to quantum 
dot itself (rather than the point contact) to a fictitious 
lead with many weakly coupled channels 
As shown in appendix O this model results in an effec- 
tive action S similar to H3.17|l . but with a sub- unitary 
scattering matrix S. If relaxation is strong, S is small, 
and expanding in the action is justified. For this situa- 
tion, we recover the interaction corrections found in Refs. 



B. Large number of channels 

We have not been able to do a full calculation to all 
orders in S. However, by organizing each order in the 
effective action to the (formal) power of TV it carries, we 
have been able to calculate the current to all orders in S 
while expanding in 1/A'^. The reason that an expansion in 
1 /N is possible is that the correlator (|3.35(l of interacting 
fermions admits a systematic expansion around the non- 
interacting correlator, by expanding 

/i/^ = l + lln/ + .... (4.11) 
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Counting the power ofl/N is done keeping the RC time 
tt/E^N and the dwell time Td constant, so that the only 
factors arise from the expansion H4.11(l and from ex- 
plicit summations over the channel indices. Oscillating 
contributions with an explicit dependence on the dimen- 
sionless gate voltage J\f, involve A^th order in perturba- 



tion theory and cannot be considered in this method. 

Non-oscillating contributions to the current occur to 
even order in the effective action only. The non- 
oscillating current contribution of order 2n in the action 
reads 



I J, nit) 



i2n) 



.(-z)^"(5"/,)o 



(-2w_f)^" ^ ... ^ f dti . . . dtndsi . . . dsn f dri . . . dvndai . . . danVk^ 

X W/i • • ■ VlNSki,k{itl,ti - Ti) . . . Sk„.k'^{tn, tn ~ Tn) (S'''' );'^ ^ (si - (Ti, Si) . . . (S'^)i;,i„ (Sn - Cr„, S„ 
X ('?'cV^iiL(^l)V'fc;i?(^l - n)l^i\fl(si - Cri)^A;^i(si) . . . 



(4.12) 



As before, the time differences denoted by Greek sym- 
bols do not affect the contour ordering for the fermion 
operators. 

The average of the product of 4?! fermion operators 
and one current operator in Eq. H4.12|l admits a stan- 
dard diagrammatic representation if the factors 1 /N In / 
in the expansion 14.11|1 of the general average H3.35|l are 
considered 'interaction lines'. Without 'interaction lines', 
the average is simply given by Wick's theorem: Denoting 
each single-fermion Green function by a solid line. Wick's 
theorem expresses the average as a sum over partitions 
in 'bubbles' of alternating single-fermion Green functions 
for right moving and left moving fermions. One bubble 
contains the current vertex i. In our formal expansion in 
1/A^, each bubble contributes a factor N, as it contains a 
trace over the channel index. Each 'interaction line', on 
the other hand, contributes a factor 1/N. Only diagrams 
in which all 'bubbles' are connected by interaction lines 
contribute to the current. This means that, to leading or- 
der in 1/A^, the relevant diagrams have a 'tree' structure, 
see Fig. |5l Every 'bubble' has precisely one connection 
to the current vertex via intermediate bubbles and in- 
teraction lines. Subleading contributions to the current 
contain 'loops'. 

To leading order in we only have to examine the 
'tree' diagrams, plus the current contribution that is ze- 
roth order in the action, see Eq. (|4.1|) . Combining these, 
we find 

I jit) = ^ / dti^ijkit -ti)[nLkiti) - ^J-Rkitl)] 



fc=i 



+ e / dhtrtjit-hyj^iti), 



(4.13) 



I(t) 



1 + 



(a) 




FIG. 5: (a) Diagrammatic representation of leading-in-l/A^ 
contribution to the current. The current vertex i is denoted 
with a filled circle, the 'bubble' of single-fermion Green func- 
tion with a solid loop, and the tree structure with a hatched 
triangle. The defining equation for the tree structure is shown 
in (b) , where the dotted arrow represents the interaction line. 



reflection matrix r. In Eq. H4.13|) . we abbreviated the 
contribution of the bubble by 



1 

J-it) = 



Gfit, t)- j dTdaSit,t~T) 
G^it-T,t-<T)S^t-cr,t)\ (4.14) 



where 



G^it, s) = G^it, s)e"^('^(*)-^(^)). (4.15) 
The phase 0(i) satisfies the self-consistent equation 

^(t) = 2 / dtiKo(ii -i)trr2j-(ti). (4.16) 



where we used Eq. (|3.14|l to express v in terms of the (|4.14|) ]. 



The quantity j_ is nothing but the Landauer expression 
for the difference of the current flowing into the leads 
coming from the quantum dot [second term in Eq. 14.14|l ] 
and coming from the fictitious reservoir [first term in Eq. 
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the dot potential is 



+ 



t N 



(4.17) 



FIG. 6: Diagrammatic representation of the eflective interac- 
tion K, (thick arrow) in terms of the bare interaction line In / 
(thin arrow) and a 'bubble' of single-fermion Green functions 
(solid loop). 



The overall phase factor that appears in G^^(t, s) can 
be written in terms of the electrochemical potential /Zc 



(t) + (j>{t) - (f>^R{s) + (f>{s) 



dt'^icm{t')- 



Here ^cm — l^ra + eVd [the plus sign arises because the 
interaction is in the leads, not the dot, c/. Eq. H3.21|l ] and 



where C is the dot's geometrical capacitance and we used 
the specific form of Ko(i' ^ t) and ijk{T). Note that in 
the limit rj 1 the (t) drop out of (|4.13l) . Equation 
(|4.17l) is just the same as (|l.lf) . written in terms of the 
currents. Hence, to leading order in 1/A'^, our formalism 
recovers the self-consistent theory for the effect of the 
Coulomb interaction on transport through quantum dots. 

To subleading order in 1/iV, one has to solve diagrams 
with one loop. The 'interaction lines' in these diagrams 
are replaced by an effective interaction, which is an RPA- 
like series, see Fig. |S| Denoting the effective interaction 
by K(t, s), one finds that K{t, s) satisfies the self-consistent 
equation 



K{t,s) — KQ{t,s)~2vp J dtidTidai{K(){t,ti ~ Ti) ~ KQ{t,ti ~ ai)) 

X trr''S{h,h - Ti)Gl{h - Ti,h - ai)SHh - cTi,h)K{ti, s). 



(4.18) 



Note that K{t, s) — ii t > s and that K{t, s) is real. 

The one-loop interaction correction to the current can be represented by four diagrams, see Fig. \7\ The first and 
second diagram in Fig. [T] represent a correction to j_. 



Gfit, t)~ J drdaSit, t ~ T)G\(t -T,t-(T)S\t- ct, t) 
u|^Im j dsidTidaiK{si,t — Ti)S{t,t — Ti)G^{t — Ti, si — ai)S'^ (si — ai, si) 
(1 - r^)Gf{si,t) +r^ I dT2da2S{si,si - T2)G1(si - T2, t - a2)SHt - a2,t) 



(4.19) 



whereas the third and fourth diagram represent an additional renormalization of the distribution function for right 
moving fermions, which can be represented by a change in the relation between to Gr and Gr, 



G'^{t,s) = Gl(i,s)e-^('^W-^(^))+"(*.^). 



(4.20) 



For the calculation of a(t,s), the rcgularization of the fcrmion operators at the impurity site is important, which 
makes the actual calculation rather cumbersome. The final result reads 



a{t, s) 



1 










(1 


X 


(1 



dtidt2{K{ti -t)- K{h - s)){K{t2 -t)- K{t2 " s))tr \ Gf {h , t2)Gf {t2 , h) 

(T2)S\t2 ~ (J2,h) 



(4.21) 



Equations H4.ia|l . H4.16|l . H4.18|l . H4.19|l . H4.20|l . and (14.21(1 give the solution of the problem up to sub-leading order 
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c 



c 




FIG. 7: Diagrammatic representation of subleading-in-l/A'^ 
corrections to the current 



in l/N. 

Note that the leading corrections to the theory are 
in the form of a modified expression for the current j_; 
Equation (|4.16|) and, hence, the relation (|4.17|) between 
Ij{t) and the self-consistent potential Vd remains un- 
changed. 

The interaction correction to j_ , contained in the sec- 
ond and third line of Eq. H4.19|l . consists of a term pro- 
portional to and a term proportional to 1 — . When 
expanding in the action and truncating the expansion at 
second order, as is done in Sec. II V Al and in Refs. i2i and 
1^ one recovers the first term only, with prefactor unity 
instead of 1 — r^. However, upon setting — > 1, which 
is the relevant limit for a coherent quantum dot with a 
unitary scattering matrix, the prefactor of that first term 
vanishes, and the interaction correction to j_ is given by 
the second term instead. It is this replacement of the 
first term by the second term in the interaction correc- 
tion that is the essential difference between the theories 
of Refs. H and on the one hand and that of this work 
and Ref. l2J on the other hand. 



where A was defined in Eq. p. 4(1 above. We will be using 
the Fourier transform of the scattering matrix and the 
Green functions. 



S{e) 
S\e) 



dtS\t)e''\ 

= —Smn{2fR,Lis)-l), 
VF 

where / is the distribution function. 



(5.2) 
(5.3) 

(5.4) 



A. Fully coherent dot 

We first consider the case = 1, corresponding to a 
fully coherent dot. 

We first consider the function (j>. By Eq. H4.16|l . 4> 
depends on trj_ only. For the time independent case, 
one has 



tr j- 



(5.5) 



Calculating the derivative dcp/dt from Eq. (|4.16() one then 
finds 



V. STEADY STATE TRANSPORT 

In this section we analyze the results of the previous 
section for the case of steady state transport through the 
quantum dot. We assume that the chemical potentials 
in the reservoirs and the gate voltages defining the shape 
of the dot are time independent. This implies that the 
scattering matrix S(t, t — t) is independent of t and the 
Green functions GL{t,t — T) and G}i(t, t — r) are indepen- 
dent of t. Current conservation then implies that there is 
no net current into the dot, /j = - see Eq. H3.27|l . 
This allows us to focus our calculation on the weighted 
difference of currents in the left and right leads, 

N 

I = J2A,Ij, (5.1) 



9l _ _N f dKo{t-ti) d} 

dt ~ TT J ' dt dt 



which implies that the derivative d(j)/dt is left undeter- 
mined by the self-consistency equation (|4.16() : any time- 
independent value for d(j>/dt is a solution of the self- 
consistency equation. This is no problem for steady state 
transport; in the previous section we have seen that (jj is 
determined once the time-dependence of bias and gate 
voltages is taken into account. In this section, we will set 
= from now on. 

The effective interaction function k can be solved from 
Eq. H4.18|l using Fourier transform, 
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k{uj) ~ K,a{to) 



^.aiiu) - - / detvSie - Lo){fR{e - lu) - fR{e))S\e) 



Here kq{u!) is the Fourier transform of kq{t), 



ko{llj) = / dTKo{T)e 



(5.7) 



This effective interaction is the same as that obtained by Golubev and Zaikini^ We further need to calculate the 
hmction a{t — s), see Eq. ()4.2U|I . Fourier transforming Eq. (|4.21|) one has 



(5.8) 



v{t) = J duj{l - cos{ljt)Muj)\^ 

jNcoth{uj/2T)-J2 J detrS{e±uj)fR{e±Lj)S^{e±uj)S{e){l- fR{e))S\e) 



Finally, with these results the current through the dot is 

/ = --^trA/ij^ + -^ / detYAS{e)fR{e)S\e) 



Im / dujdeK{uj)tYk[S{e)fR{e + uj)S'^{e + uj)S{e){l- fR{e))S^{e) 



S{e){l - fRie + Lo))S\e + to)S{e)fR{e)S\e)] 

du;de\K{Lu)\hYAS{e)[2fR{e) - /^(e - uj) - fR{e + co)]S^{e) 



^TVcoth^-^ J de'trS{e'±u;)fR{e'±uj)S\e'±io)S{e){l^fR{e'))SHe') 



In order to find the linear conductance, we set fiR = —AeV, fiL = 0, expand in V, and find / = GV, with 



(5.9) 



G = ^T-rtrA^- — / de 



Im / deduj 



-^iP-\trKS{e)KSHe 
OS I 



(5.10) 



27r2n 



de 



K{Lo)tv [KS{e)KS\e + lo) - KS{e)S\e + uj)S{e)kSHe)\ (1 - 2f{e + lj)). 



The same equation was derived previously by Golubev and Zaikin using a different method^ 

Equation H5.9|) gives the current for one particular realization of the scattering matrix S{e). In order to find the 
average and variance of the current or the conductance for an ensemble of quantum dots, we need to average 5'(e) 
over the appropriate ensemble of scattering matrices. Details on the averaging procedure can be found in Appendix 
IbI For large N , the average of a product of traces may be calculated as the product of the averages. This means that 
the interaction kernel k{uj) may be averaged separately, with the result 



1 + ILOT^ 



(5.11) 



In Eq. (|5.11(l . Td is the mean dwell time for electrons entering the dot. In terms of the dot's mean level spacing A 
and the total conductance <? = 5i + 52 of the point contacts, one has 



27r 
5A- 



(5.12) 



(If levels are spin degenerate, A is half the mean spacing between spin-degenerate levels.) Performing the ensemble 
average for the current, one then finds, to leading and sub-leading order in N and in the absence of time-reversal 
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symmetry, 



" -^trA^ifl + ^y de[gtrArJ(e)rt+trA(l-r,rt)tr/(e)(l-rtrJ] 

/(iujT )^ 
deduiniu})— — (5.13) 

X [gixKrcf{e + Lu)rl{l - r,rj) tr(l - /(£))(1 - rlr^) 
+ 5tr Ar,(l - /(e + w))rt(l - r^rj) tr /(e)(l - r^r,) 

- tr A(l - r,rt) tr f{e + ^) (1 - rjr Jrjr, tr (1 - /(e))(l - rjr J 

- trA(l - r,rt) tr(l - f{e + u)){l - r^r^rtr, tr /(e) (1 - rjrj] , 

From Eq. H5.13|l we conclude that there is no interaction correction to the average current for ideal point contacts. 
Setting = —eVK and expanding Eq. \h.VA\ to first order in the bias voltage V, we find the ensemble average of 
the linear dc conductance, 

(G> = (l - Im / do. (g2.si+ffi.2)r,[(^/r)~sinh(a;/r)] 1 

2Tih g \ J 2g2(l + icjrd)[l + (7r/i;e5Td)(l + iwTd)]sinh2(cj/2r)J 



Performing the frequency integral in the limit K/t^ <^ 
T < E^g one findsSi 



(G) 



9i92 



2'Kh g 



1 



2(si52 



S2gi) £^c5e 



i+c 



r 27r2r 

(5.15) 

which is the same result as for a quantum dot without 
coherent scattering from inside the dot, see Eq. (|2.8|l . In 
the low temperature limit T <gC S/ra <C Ecg one finds a 
saturation of the interaction correctioni2i 



(G) = 



2iTh g 



2(52Si 



9' 



(5.16) 



For the quantum interference corrections we have per- 
formed calculations for the linear conductance with ideal 
point contacts only. Weak localization is a small differ- 
ence 5G of the ensemble averaged conductance with and 
without magnetic field. We find that the interaction cor- 
rection to weak localization is zero, so that 5G is given 
by the non-interacting theoryj^SiS 



5G = - 



Vse^NiN2 
2nhN^ 



(5.17) 



where i>s denotes the spin degeneracy. 

Closed form expressions for the interaction correction 
to the conductance fluctuations could be obtained for 
the limiting cases T < fi/rd < NE^ and fi/rd < T < 
NEc, for which we find, in the absence of time-reversal 
symmetry 



var G = 



2TrhN^ 



E,N^ 



1 - ^ {nTTdf 



(5.18) 



and 



varG 



errd 



respectively. For intermediate temperatures, the final re- 
sult still contains two energy integrations, which are eas- 
ily integrated numerically. Figure|Slshows the interaction 
correction to the variance as a function of temperature 
for the value n/EcNrd — 0.1. In the presence of time- 
reversal symmetry, the variance of the conductance a fac- 
tor two higher than in the absence of time-reversal sym- 
metry. We conclude that, at zero temperature, there is no 
interaction correction to the conductance fluctuations, to 
leading order in 1/N. At small but finite temperatures, 
interactions lead to a small decrease of the conductance 
fluctuations. 

It is legitimate to ask how one can distinguish between 
the interaction correction to weak localization and the 
conductance fluctuations and the second order interac- 
tion corrections to the classical [i.e., of order N] conduc- 
tance or the second order quantum interference correc- 
tions in the non-interacting theory, as both are of order 
1/A^. With respect to the comparison to the second or- 
der interaction correction to the classical conductance, an 
answer is readily given: this interaction correction does 
not depend on the magnetic field and has no mesoscopic 
fluctuations as a function of, e.g., Fermi energy or mag- 
netic field. A formal method to distinguish the interac- 
tion corrections discussed above from second-order quan- 
tum interference corrections follows after introducing the 
number of 'orbital channels' N° = N/vg- Then one finds 
that the interaction correction is a l/N correction, while 
the second-order quantum interference correction is of 
order N/{N°)'^. For the conductance fluctuations in the 
absence of time-reversal symmetry this problem does not 
arise, as there are no second order quantum interference 
corrections in this case42*^Si^ 
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27tT/NA 



FIG. 8: The interaction correction to the variance of 
the conductance, made dimensionless by writing G = 
{vse'^/2-Kh){NiN2/N^)g. The solid curve shows the non- 
interacting contribution to the conductance fluctuations; the 
dashed curve shows the interaction correction for -k/ E^Nt^ — 
0.1, multiplied by to make it fit on the same scale as the 
non-interacting contribution. 



theory of Refs. HEIl For the linear conductance one then 
finds 



2 r 



G 



e 
2^ 



(5.21) 



if Iii{Ec/ttT) < N/2 but T < E^N and 



G = 



2TTh 



N1N2 

N 
E^Ne^ 

TT^T 



27ri/2 I N \N 2 



A ? t " 

tr A r„rJ. cos — 

C C 



(5.22) 



if ln(£'c/7rr) ^ N/2 (but temperature larger than a suit- 
able lower limitSiiiS) . 

For large N, one recovers the same expressions as Gol- 
ubev and Zaikin, see Eq. (|2.8|l and Ref. 123, starting from 
the results of Sec. HVbI 



The absence of an interaction correction to weak local- 
ization, and to the conductance fluctuations at zero tem- 
perature, is perhaps one of the most striking predictions 
of our theory. In Section rVIII we will give an argument to 
explain why such a renormalization leaves these ensem- 
ble averaged quantities unchanged, using the observation 
that the form of the conductance correction derives from 
the renormalized scattering matrix Eq. (|3.6|l . 



B. Incoherent dot 

The opposite limit of a quantum dot without coherent 
reflection from inside the dot can be obtained from the 
general formahsm of Sec. HVl bv setting S{t,t — T) — > S{t) 
and interpreting r as the reflection matrix r^, of the point 
contacts. In this case, the distribution function of the left 
moving fermions represents the distribution function of 
electrons leaving the quantum dot. Self-consistent deter- 
mination of the corresponding chemical potential /ii is 
problematic, however, because that would require knowl- 
edge of the average dwell time inside the dot, which is not 
contained in this model. 

For small Tc, the results of Sec. IIV Al can be used. In 
this case one flnds 



ZTTft. 

1 f°° 

\- -eivp / dTsin(Ko(-T)) 

2 Jo 



(5.20) 

l/N 



f{r)f{-r) 

X t,Kr,rl{Gl{T)Gl{-r) GI{t)GI{-t)). 

For N ^ 2 there is an additional oscillating interaction 
correction to the current. Equation (|5.20|) is the non- 
equilibrium generalization of the original linear response 



VI. TIME-DEPENDENT TRANSPORT: 
ADIABATIC APPROXIMATION AND LINEAR 
RESPONSE 

In this section we will consider the case of a fully co- 
herent dot only. We will limit ourselves to a theory of 
adiabatic and linear response transport — the response 
to either a slowly varying internal potential or to a small 
bias voltage, but not both — , and consider the average 
and variance of the current only. 

We will present our final results in terms of the Fourier 
transform of the scattering matrix, 

S{t;e) = J dTe"'^Sit + T/2,t-T/2), 
S''{t;e) ^ I dTe'^-^SHt + T/2,t~T/2). (6.1) 



In the adiabatic approximation, the Fourier transform 
S{t;e) satisfies the unitarity condition 



il-{i/2)Dt^,)S{t;e)S\t;e)^l, 



(6.2) 



where Dt^.AB = {dA/dt){dB/de) - {dA/de)[dB/dt). 
To leading order in the adiabatic approximation, the 
Fourier transform S{t\ e) is equal to the 'frozen' scatter- 
ing matrix, taken by fixing the internal potentials at the 
value they have at time t. However, the 'frozen' scatter- 
ing matrix is unitary, whereas S{t] e) is not. Since our 
final expression for the current will be manifestly of first 
order in combined linear response and adiabatic approxi- 
mation, we can neglect the difference between S'(i; e) and 
the 'frozen' scattering matrix in our final expressions. 

We will set the chemical potential hl of the electrons 
in the fictitious reservoir equal to zero for all times. In 
the adiabatic approximation and for linear transport, the 
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Keldysh Green function Gj^jj.(t, s) depends on the time 
difference t — s only, m — 1, . . . , N, and reads 



Gl„j,{t + T/2,t-T/2) 



T 

VF 



X P- 



(6.3) 



-i(p™H(t)+eVd(t))r 

sinh(7rrT) 



where Vd is the dot potential, see Eq. (|4.17() . 

Calculating the interaction function n and the inter- 
action correction a in the adiabatic limit is problem- 
atic. The reason is that both quantities involve long 
time scales during which one cannot assume that the 
derivative of the scattering matrix is constant (as one 
does in the adiabatic approximation). However, since n 
and the interaction correction a contain traces, they are 
self-averaging, and we can replace them by their ensem- 
ble averages. Taking the ensemble average in the defin- 
ing equation (|4.18l) for K{t, s), using the known probabil- 
ity distributions of time-dependent scattering matrices, 



the equation for nit, s) becomes 



EqN f /"'^ 

K{t,s) = Ko{t — s) -\ / dt2 / dt'2K(t2,s) 

Jt Jt 



(6.4) 



This equation has no reference to the distribution func- 
tions in the lead or the time dependence of the scattering 
matrix. Hence we find that the solution is of the equi- 
librium form K{t — s), the Fourier transform of which 
is given in Eq. H5.11|l . For the interaction correction 
a{t + T/2,t — t/2), we note that it is an even function of 
T. This implies its ensemble average must vanish if we 
consider the adiabatic approximation and linear response 
only. 

We are now in a position to calculate the self-consistent 
potential Vd and the current Ij . The self-consistent equa- 
tion for the potential Vd reads 



dt n J \ de J I dt de dt 

+ iRe j du^K{oj){2f{e + uj- t,,)-\)(^^^ + ^^[S{t-e)S\e + oj)]Y (6.5) 
where the distribution function / is evaluated with respect to the electrochemical potential 

k 

The first terms give the mean- field contribution to Vd. The last term gives the interaction correction. Upon taking 
the ensemble average, which is appropriate if Vd enters as a variable in the expression for the current, the last term 
vanishes. 

For the current we find 



CdV, e [ ( df{e-j^\ ( 1 
N^^2^r'\ d^)'' V^''N 



X {j{t;e) + -lT[v I dujK{uj) {2J{e + uj - 11,) - 1) S{t-e)S\t-e + uo){j{t-e + u) - j{t-e)) \ , (6.7) 



TT 



where we abbreviated 

j{t-e) = ^j,^S{t;s)(^nn-t^y^^t-I^^SHt;e). (6.8) 

In the expressions for Ij{t) and dVd{t)/dt one may replace the Fourier transformed scattering matrix S{t;e) by 
the 'frozen' scattering matrix. Then a time derivative d/dt is replaced by {dX/dt){d/dX), where X is a (set of) 
parameter(s) that characterizes the potentials in the dot. 

In Eq. H6.7|l the terms proportional to /i^ give the dc conductance of the dot, cf. Eq. H5.10|l . The first line of Eq. 
(|6.7fl is nothing but the Landauer formula with a self-consistent potential in the dot; the last lines give the correction 
from interactions beyond the mean-field level. The self-consistent potential Vd was omitted from Eq. H5.10|l because 
it cannot be determined from the steady state formulation of the problem only. The terms proportional to the time 
derivative give the emissivity —dQ/dX of the quantum dot. Again, the first two lines of Eq. (|6.7|l correspond to the 
emissivity according to the self-consistent (Hartree) theory of Biittiker and coworkers^ whereas the remaining lines of 
that equation give the correction from interactions beyond the mean-field level. The governing equation H6.5|) for the 
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self-consistent potential Vd (t) is the same as that of the Hartree theory of Ref. |23 It depends both on (changes in) 
the dot potential and the chemical potential in the lead. After taking the ensemble average there are no interaction 
corrections to the Hartree theory for Vd- 

Although the above expression for the current was derived in linear response and in the adiabatic approximation, 
it can be used to calculate the electrochemical capacitance Cc of the quantum dot, the derivative of the charge on the 
dot to a uniform shift if chemical potential /i^j. This is possible because, for the Hamiltonian Hq, slowly increasing the 
chemical potential of all channels simultaneously corresponds to a time-independent small shift of the electrochemical 
potential, see Eq. H3.32|l . Hence, the time derivatives of the phases carried by the bare Green functions Gr and Gl 
differ only by a small amount, and linear response is meaningful. Considering the rate of change of the charge on the 
dot for the case in which there is no explicit time dependence of the dot scattering matrix, one finds from Eq. (|6.7|l 



dt 



dt ' 



(6.9) 



where dVd/dt can be calculated from Eq, 
taneously by an amount —eV{t), we findS& 



Upon changing all chemical potentials for incoming electrons simul- 



- dV 



= [G-^ + {e^dn/de)-^] \ 



(6.10) 



where dn/de is the density of states of the dot 
9/(£- Mc) 



dn 1 
d^ = 2^''^' 



de 



tr 



iSis)^^l^ + -Re / dL0Kiu)i2fis + u;-^l,)- l)-^Sie)S\e + lo) 
oe IT I oe 



(6.11) 



Here, the first term is the Hartree expression for the den- 
sity of states of an open quantum dotj^Sj^l the second 
term is the interaction correction. 

The ensemble average of the interaction correction to 
the density of states, and hence of the interaction correc- 
tion to the electrochemical capacitance, is zero. It has, 
however, mesoscopic fluctuations. Calculating the inter- 
action correction to the variance of the density of states, 
we find, for zero temperature. 



var [dn/de) 



27r^ 
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from which one concludes, for h/rd ^ EcN , 



(6.12) 




27rT/NA 



var Cc 



2EtN^Tl 
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(6.13) 



Equations 1)6.12(1 and (|6.13|) are valid in the absence 
of time-reversal symmetry. With time-reversal symme- 
try, the variance is a factor two larger. The first term 
is the variance of the mean-field contribution to the 



capacitance 



36.58 



the second term is the interaction cor- 



rection. Figure shows the temperature dependence of 
the variance of the electrochemical capacitance. 

If the quantum dot is operated as a quantum pump, 
two shape-defining parameters Xi and X2 are varied pe- 
riodically in time, whereas the bias voltages are kept zero. 
We consider a harmonic time dependence for the param- 



FIG. 9: The variance of the mean-field contribution 
to the capacitance, made dimensionless by writing Cc = 
2-Ke'^c/{E^N'^Td) (solid curve) and the interaction correction 
to the variance of the electrochemical capacitance (dashed), 
as a function of temperature. The interaction correction to 
the variance has been multiplied by A'^ in order to fit on the 
same vertical scale. 



eters Xi and X2 



Xi{t) = SXisin{ut), 
X2{t) = 5X2fi\n{uJt + 



). (6.14) 
is defined as the 



The charge pumped in one cycle Qpump 
integral of the current in one of the point contacts. Inte- 
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grating over one cycle, one then finds, for small pumping 
amplitudes'- 



?pump = {716X16X2 sin(f>)^ 

dx'i dX2 'dX^ dXi 



(6.15) 



[For large amplitudes 6X1 and 6X2 one has to average 
Eq. H6.15|l over the area enclosed in the two-dimensional 
space spanned by the parameters Xi and X2.] 

Since the total current is conserved, we may calcu- 
late Qpump as the integral of I{t) = J2j O'^'^r 



one period, i.e., {N2/N) times the integral over the cur- 
rent through the first point contact minus {Ni/N) times 
the current in the second point contact. The ensemble 
average (Qpump) is zero for symmetry reasons. When 
calculating the variance of the pumped charge, we note 
that, up to corrections of order l/N'^, one may replace 
the average of a product of traces by the product of the 
averages, with the exception of the two traces that con- 
tain the matrix A. Their average is zero, and one needs 
to consider the average of the product. This means that 
all terms proportional to dV^/dt can be dropped, and 
one may set 



dQ 



pump 

dx 



e 
1 



de 



dm 



de 

Re / dujK(uj){2f{e 



tr <^ iAS{e) 



- LO 



dX 
l)trA5'(e) 



dS\e + uj) 



dX 



S\e + uj)S{e) 



dX 



(6.16) 



r 



Using the results of appendix IbI to perform the ensemble 
average, we find that, at zero temperature and for small 
pumping amplitudes, the variance of the pumped charge 
is given by 

vargpu,„p = ^^^^^ [l + ^) {6X^6X2Bin4>fiQ.n) 

irrespective of the presence or absence of a time-reversal 
symmetry breaking magnetic field. Here we have used 
the standard convention (within random matrix theory) 
to relate the effect of the shape-defining parameters Xi 
and X2 to a parametric motion of the scattering matrix. 
The first term in Eq. H6.17|l is the leading non-interacting 
contribution to the variance >2fliS£ the second term is the 
interaction correction. For temperatures T ^ h/rd we 
find 

varQpump = {-K6Xi6X2S\n(l)f 

where, again, the first term is the leading non-interacting 
contribution to the variance of the pumped currenljSSiSli 
and the second term is the interaction correction. 



VII. CONCLUSION 

In this paper we have presented a systematic theory 
to evaluate the effects of electron-electron interactions 
on quantum transport through open quantum dots, ac- 
commodating arbitrary bias and time dependence of the 
dot's scattering matrix S{t^ t'). Our result takes the form 
of an expression for the sample-specific current through 



a quantum dot in terms of S{t,t'). This expression 
has been obtained using a systematic expansion in 1/iV, 
where N is the total number of channels in the point 
contacts connecting the dot to the electron reservoirs. 
To leading order in 1/7V, our theory reproduces what 
one finds if the interactions were treated in the Hartree 
approximation. The subleading-in-l/7V corrections rep- 
resent contributions to the current through the dot that 
cannot be described by means of a self-consistent poten- 
tial. 

Upon taking the average over an ensemble of chaotic 
quantum dots, we have been able to calculate the inter- 
action corrections to weak localization, universal conduc- 
tance fluctuations, capacitance fluctuations, and to the 
pumped current in case the quantum dot is operated as 
a quantum pump. For a quantum dot with ideal point 
contacts, we found that the interaction correction to the 
weak localization correction to the conductance vanishes, 
as well as the interaction correction to the conductance 
fluctuations at zero temperature. There is a negative in- 
teraction correction to the conductance fluctuations at 
flnite temperature. There are positive interaction correc- 
tions to the capacitance fluctuations and to the fluctua- 
tions of the pumped current, but no corrections to their 
average. All interaction corrections are small as 1/A^ in 
comparison to the non- interacting {i.e., Hartree) contri- 
butions. 

The problem of 'weak Coulomb blockade' in open 
quantum dots has been addressed previously in the 
literature ^iSi^ Our flndings for the quantum interference 
corrections, as well as our expression for the sample- 
speciflc dc conductance and capacitance, differ from 
those of Refs. 0S To explain the difference, we have 
shown that the formal expansion in the scattering ma- 
trix, used in Refs. 00, cannot be used to describe a co- 
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herent quantum dot. We believe that the expressions 
presented here, which were obtained using a systematic 
expansion in are correct. Our sample-specific ex- 

pression agrees with that obtained in Ref. 24, although 
Ref. 1^ does not consider the quantum corrections to the 
conductance. 

The absence of an interaction correction to weak lo- 
calization and to the conductance fluctuations at zero 
temperature can be understood using a simple argument: 
Since the charging interaction cannot lead to dephasing^ 
— it corresponds to a time-dependent uniform shift of 
the dot's potential — , its only effect is a renormalization 
S{e) S'{e) of the scattering matrix S{e) of the dot. 
To leading order in perturbation theory, such a renormal- 
ization was calculated explicitly in Sec. IIII Bl and in Ref. 
124 . Like the probability distribution functional P[S'(£)] of 
the 'bare' energy-dependent scattering matrix S{e), the 
probability distribution functional P[S'{e)] of the renor- 
malized scattering matrix S'{e) is invariant under left 
and right multiplication of 5"(e) with energy- independent 
unitary matrices U and y^SliSiSMi 

P[S'{e)]=P[US'{e)V]. (7.1) 

In the presence of time-reversal symmetry, U = V"^ . The 
statistical invariance H7.1|l . together with the same invari- 
ance for P[S{e)], implies that the distributions of both 
S{e) and S"(e) for a fixed energy e are those of the circu- 
lar ensemble from random matrix theory. "^"'^^ Since the 
weak localization correction to the conductance and the 
conductance fluctuations at zero temperature do not de- 
pend on statistical correlations between scattering matri- 
ces S{e) at different energies, we conclude that these are 
given by the circular ensemble averageSf^SiS and, hence, 
that there is no interaction correction. There is an in- 
teraction correction, however, to the conductance fluc- 
tuations at a finite temperature, which do involve such 
correlations. The above argument does not apply to a 
quantum dot with nonideal leads, for which the invari- 
ance property (|7.1|l does not hold, or to a quantum dot 
with relaxation, for which the scattering matrix is sub- 
unitary and the invariance property (|7.1(l is not sufficient 
to fix the distribution. 

How reliable is the expansion in and how signif- 
icant are the interaction corrections, as they are a fac- 
tor 1/N smaller than the Hartree contributions? Re- 
calling that N is the total number of channels in the 
point contacts connecting the dot and the electron reser- 
voirs, not counting spin degeneracy, one has = 4 for 
a dot with two spin-degenerate single-mode point con- 
tacts. Although = 4 cannot be classified as 'large', 
we believe that for = 4 the expansion should 
give a reliable estimate of the interaction correction, 
whereas, on the other hand, the interaction corrections 
should still be significant (if they are nonzero). Mea- 
surements of the conductance distribution in single- mode 
quantum dots have shown excellent agreement with the 
(non-interacting) random matrix theory down to the low- 
est temperaturesJS; This experimental observation is in 



agreement with our theory since, to leading order in 1 /N 
and at zero temperature, we find no interaction correc- 
tion to the average and variance of the conductance. The 
argument following Eq. 1)7. 1|) is extremely suggestive and 
may hint that the agreement with random matrix the- 
ory is not restricted to the approximations of this paper, 
though we caution that higher orders in the interaction 
are presumably not simply described by an effective S'{e) 
due to inelastic processes. As long as the Fermi Liquid 
picture holds, however, one may expect this description 
to be restored at zero temperature. 

The interaction correction we calculated here is the re- 
sult of the capacitive interaction only. For large N, it is 
a factor ~ 1/iV smaller than the leading non- interacting 
(Hartree) contribution to the transport properties and 
their quantum interference corrections. We have not con- 
sidered other contributions to the electron-electron in- 
teraction. These were omitted from the Hamiltonian of 
the quantum dot because they are much smaller than 
the capacitive interaction if the dimensionless conduc- 
tance of the dot is large&SiS^ The interactions we omit- 
ted will give additional corrections to the dot's transport 
properties, e.g., because they lead to dephasing. To the 
best of our knowledge, a microscopic theory of the effect 
these residual interactions on the transport properties of 
a chaotic quantum dot still has to be developed. 

In addition to the interaction corrections studied here, 
which are not periodic in the dimensionless gate volt- 
age J\f, there exist periodic-in-A/" interaction corrections 
to the dot's capacitance and conduct ance&L^ These oc- 
cur to Nth order in perturbation theory, which rules out 
a large- Af approach. Periodic-in-gate-voltage interaction 
corrections have been studied using an expansion in the 
scattering matrixA^ Although we have shown that such 
an expansion is not justified for aperiodic interaction cor- 
rections, we must leave the question of its validity for the 
periodic-in-A/" corrections unanswered. 

Before closing, we would like to make two remarks on 
the relation of our work to that of Golubev and Zaikini24 
First, Golubev and Zaikin consider both the case of a 
'voltage biased' and a 'current biased' quantum dot. The 
latter case corresponds to a quantum dot in series with 
a large ohmic resistor. In Ref. it is shown that the 
interaction corrections in the two cases are qualitatively 
different. Our calculation, in which the chemical poten- 
tials of the electrons coming in through the point contacts 
are kept fixed, corresponds to the 'voltage biased' case. 
Second, Golubev and Zaikin argue that the sole effect of 
interactions as a renormalization of the reflection matrix 
of the point contacts Tc only. While this is true for an 
incoherent dot,^ii^2iSLS our calculation shows that this 
is not correct for a coherent dot: When quantum inter- 
ference corrections are taken into account, the renormal- 
ization of the scattering matrix is more complicated, and 
involves more than one energy. This is manifested, e.g., 
in the existence of an interaction correction to the con- 
ductance fluctuations at flnite temperature. 
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APPENDIX A: BOSONIZATION 

In this section we calculate the contour-ordered corre- 
lation functions of the fermion fields in the Hamiltonian 
Hq of Eq. (|3.18() . The bosonization method allows for an 
exact calculation of the fermion correlation functions for 
the Hamiltonian Hq, despite the fact that Ho is quadratic 
in the fermion operatorsiSS Hq is rewritten in terms bo- 
son fields (pjL{x) and (pjR^x). In the bosonized language 
Hq is diagonalized easily, and the fermion correlators can 
be found from the boson correlation functions. 

The boson fields if>jL{x) and Lpjji{x) are related to the 
fermion fields ipj^lx) and ^jji{x) as 



There are no boundary terms at a; = — cxd because the 
interaction extends over the entire lead. The reference 
electron number A'lcf has been absorbed in the definition 
of normal-ordering of the fermion operators. Then, once 
expressed in terms of the boson fields, the Hamiltonian 
Hq becomes quadratic, 



Ho 



VF 

in 



N 

E 

47r2 



dx 



\ dx 



dx 



N 



^(^,L(0) + ^,K(0) + 27rAA) 



(A8) 



e-'^^\ j = l,...,7V, (Al) 



They obey the commutation rules 



.N. 



(A2) 



The Hamiltonian is supplemented with the condition 
that the chemical potential of left-moving fermions mov- 
ing towards the dot-lead interface at x = is equal 
to iJLjL(t) /2t:vp and that the density of right-moving 
fermions is equal to fijii{t)/2TTVF- We generate the elec- 
tron densities corresponding to these chemical potentials 
by the inclusion of a time-dependent forward-scattering 
potential in the Hamiltonian. This potential is located at 
a distance a from the dot-lead interface at x = 0. (The 
precise value of a is not important for the subsequent 
considerations, as long as a > 0.) In boson language it 
has the form 



[ipjLix),(pLiiy)] = -i^^S^JsigI^{x-y), 
[(PjR{x),(pm{y)] = i7rdysign(a; - ?;), 
[<fjBXx),<fiLi{y)] = -iTrSij. (A3) 

In Eq. IIA2I) fij — fi] is a. Majorana fermion, I??,, — , , • , , r ;i i • , 

„f Tx • 1 ££ J.- J. • J. 1 J. \l \^ where m,L and TOoo are the mtegrals of the chemical 

zOij, and A is a cut-on time that is taken to zero at the . / , . , 



3 = 1 



dx 



(A9) 



end of the calculation. Since the Majorana fermions do 
not enter into the Hamiltonian, their average is given by 



potentials HjL and ^^t?, respectively, 

t-\-a/vf 



{TcVjit), %(s)) = %sign^ (t - s) 



(A4) 



where the sign sign^, is defined with respect to the order- 
ing along the integration contour c. 

The (normal-ordered) densities of left and right moving 
fermions are proportional to the derivatives of the boson 
fields, 



TOjL(i) = 
irijRit) = 



t-\-a/vF 



dt'ii,B{t'), j^l,...,N. (AlO) 



:V']L(a;)V^jL(a;): 



J_d_ 

2tt dx 
]_d_ 
27r dx 



^odx), (A5) 
^,r{x). (A6) 



Hence, the number of electrons on the quantum dot is 
simply expressed as the sum of left-moving and right- 
moving boson fields at a; = 0, 



We assume that the system is in equilibrium, fJ-jpit) — 
fJ-jR(t) = 0, J = 1, . . . , iV for times t smaller than a refer- 
ence time tj-cf. 

In order to find correlation functions of fermion cre- 
ation and annihilation operators, expressions of the form 
TcC'^ need to be averaged over the boson fields, where 
a is an arbitrary linear combination of the boson fields 
ipjL and fjR, j — 1,...,N, at different times. Since 
the Hamiltonian is quadratic in the boson fields, such an 
average can be calculated as 



(exp(a))=exp<{(a) + i[(a2)_(a)2 



dot 



1 ^ 



(All) 



2tt 



which means that only contour-ordered correlators of up 
to two boson fields at different times are needed. Hence, 
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it remains to calculate the correlation functions of the 
boson fields. 

For this calculation, it is necessary that we transform 
the boson fields such that the charge mode is separated 
from modes that do not charge the quantum dot. The 
corresponding basis change reads 



3 



(A12) 



where the transformed fields carry a Greek index and the 
original fields carry a roman index. Here the o^j are real 
numbers satisfying 



where 



3 3 

(A17) 



In the transformed basis there is no charging interac- 
tion for the modes /i = 2, . . . , TV. Solving the Heisenberg 
evolution equation for the fields ^p^L and <^^_r, 



and 



Oij 



,iV. 



(A13) 



With this choice, the transformed fields (f>iL and (pm 
describe a mode that charges the quantum dot, whereas 
all other fields ip^L and 0fj,R, n = 2, . . . , N, do not involve 
charging of the dot. The transformed boson fields i^^l 
and ipfj,R satisfy the same boson commutation relations 
as the original fields ipjL and ipjR, 

[^nL{x),(pi,L{y)] = -iTr6^^sign{x-y), 
l'ft^Rix),'P,,R{y)] = in6ij,^sign{x-y), 
['fiij,R{x),iPi.L{y)] = -iTr6^^. (A14) 

In terms of the transformed fields, the boson Hamiltonian 
Ho and the time-dependent perturbation H' read 



H' = 
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dx 



JV 



dip^iRix) 



dx 



(A16) 



d(pi^L{x,t) d(p^L{x,t) 
= VFmf,L{t)5{x-a)+VF^^ , 



diPi^Rix, t) 



VF'mi^R{t)6{x + a)-VF 



dLp^L{x,t) 



dt T / . dx ' 

we obtain an expression for if>{x, t) in terms of boson fields 

at the reference time t^ef, 



V'ixL{x,t) = ip^_,Lix + VF{t-trcl),trcl) 

+ mfj,L{t + {x- a)/vF)0{a - x), 

'Pu.R{x,t) = 'fifj,R{x-VF{t-tref),tref) 

- m^R{t -{x + a)/vF)6{a + x), (A18) 



(Recall that m^L{t) — mfiR{t) = for t < tref-) For 
the transformed mode ^ — \, the charging interaction 
is important. In this case, the Heisenberg equation of 
motion reads 



d0iL(x,t) , , ^, , d(ftiL(x,t) EcN ( ^ . , ^ 27rA/'\ 



dt dx TT \ ^yN 

d(piR{x,t) d<fiiL{x,t) EcN f ^ 27rAr'^. 

- VFmiR{t)d{x + a)-VF ^ I tpiL{0) + tpiR{0) + ] 0{-x). 



dt dx TT \ a/TV 
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and has the solution 



ipiLix,t) = ipiL{x + VF{t-t,.cf),Ucf) + miL{t+{x-a)/vF)0{a-x) 



E,N 



dt'9{-x)0{t' - t - x/wj.)e-^<^^(*-*')/'^ 



'PlR{-VF{t' - iref),ircf) " ^^^(t' - a/vp) + 01L{vF{t - t^ci),t^ci) + miL{t' - a/vp) + 

<PiR{x,t) = tpiii{x-VF{t~t,cf),t,cf)~miji{t+{x + a)/vF)0{a + x) 



2TTAf 



E,N 



dt' 



2ttN 



ifilRi-VFit' - ii.cf),ircf) - miR{t' - u/vf) + ipiL{vF{t - irof),irof) + miL{t' - u/vf) + 

(A19) 

According to Eqs. (|A18|I and (|A19p . the effect of the time-dependent potentials is to shift the boson fields by a real 
number. Such a shift affects the average of the boson fields, 



{(Pfj.L{0,t)) = m^Lit - a/vp), 



(A20) 



(^,,fl(0, t)) = -m^R{t - a/vF) - '^5^.1 + ^5^.i I dt' e-""'^^'-'"^'^ [mifl(i' - a/vp) - m,L{t' - a/vp)] 



N 



tret 



but not the connected two-point correlators. Hence, the connected two-boson correlation functions are the same as 
in equilibrium. They can be obtained by analytical continuation from the imaginary-time correlators obtained in 
Refs. 6, 8 or they can be calculated directly from the solution of the equation of motion. In the latter method one 
first calculates the spectral density A{t^t'). Since Eqs. (|A18|) and (|A19p express all boson operators in terms of the 
operators at the same reference time tref , the spectral density can be found from the equal-time commutation relations 
(|A14|) . We only need the spectral density at x = 0, 

A^L,.L(t,0 = ([93pL(0,i),(^.L(0,t')]-) 

= -i7r(5p,ysign(t - t'), 

([(^^L(0,i)^.fl(0,t')]-> 

^5,^, rdCe-WAsign(t'-t-C)" , 



J^liL,uR.it, t') 
A^R^uL{t, t') 



1 



([(^^K(0,<)(^.L(0,t')]-> 

= -1^5,, 1 - ^<5^,i / dCe-^^^^/-sign(t - C - t') 



AfiR^vR{t,t') — 



([(^^K(0,t),^.fl(0,t')]-) 
-i7r(5p^sign(t - t'). 



(A21) 



The two-boson correlation functions are then found using the fluctuation-dissipation theorem, see, e.g., Ref.rTOl The 
resulting expressions for the averages and the connected correlation functions of the original boson fields i^Lj and 
•f>Rj, j = 1, . . . ,iV are 

ft 



dt'fijLit') 



i^jRit)) ^ - J dt'iijL{t') - 
where iij- was defined in Eq. (|3.25() . and 



2ttN 
N 



/ dt'J2^jk{t - t')[fikR{t") - fikL{t" 

k=l 



(A22) 
(A23) 



{T,^^L{t')^JL{t")) = (5yln(^ 



\TTT/i 



sinh [ttT {f - t" - asign^(i' - t"))] 
{n^^L{t')^Mn) - sign,(i' - t") 1 ln/(t", t') - A 

= {T,^,R{t")ipjL{t')), 



A+\ 5i, - 



N 



B 



(A24) 
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where A and B are positive constants that are taken to 
infinity at the end of the calculation and / is given in 
Eq. H3.36|l . Calculation of the fermion correlators of Sec. 
IIII El using the rule l|Alip is straightforward now. 



APPENDIX B: ENSEMBLE AVERAGE 

In this section we present some material relevant for 
performing the averages over an ensemble of microscop- 
ically different but macroscopically equivalent quantum 
dots. 

We consider the statistical distribution of the scatter- 
ing matrix for an ensemble of chaotic quantum dots. The 
dots are placed in a magnetic field. The magnetic field 
strength is described by the dimensionless parameter a, 
which is proportional to the magnetic flux $ through the 
dot, 



e<I> 

he 



nhvpl 



-,1/2 



(Bl) 



where A is the mean spacing between spin-degenerate 
levels, L its linear size, I the elastic mean free path in- 
side the dot, and k a numerical constant of order unity. 
One has k = Air/lh for a diffusive sphere of radius L and 
K = 7r/2 for a diffusive disk of the same radiusSi. If the 
electron motion inside the dot is ballistic, but with diffu- 
sive boundary scattering, / is replaced by 5-L/8 and 7rL/4 
for the cases of the sphere and the disk, respectively. 
The shape of the dot is controlled by two gate voltages, 
which are described by the dimensionless parameters Xi 
and The normalization of the parameters Xi and 
X2 is the same as in Ref. 

Two types of scattering processes contribute to the 
scattering matrix: direct reflection at the contacts, which 
is described by the energy-independent reflection matrix 
Tc, and scattering from inside the quantum dot. The ma- 
trix Tc is a property of the contacts and is not subject to 
mesoscopic fluctuations; it is only the contribution that 



J 



involves scattering from inside the dot that has fluctua- 
tions. In order to manifestly separate the two contribu- 
tions, S is parameterized as2^ 

S = r, + {1- T,rlY'^So{l + rlSo)-\l - rlr.f'^ (B2) 

where Sq has the same statistics as the scattering matrix 
of a chaotic quantum dot ideal contacts. 

For electrons with spin (and without spin-orbit scat- 
tering), the scattering matrix takes the form 



5= 5°(^l2, 



(B3) 



where I2 is the 2x2 unit matrix (acting on the spin 
degrees of freedom) and S° is a unitary matrix of size 
N° = N/2 representing scattering of the orbital degrees 
of freedom. A similar decomposition applies to the matri- 
ces Tc and 5*0. The results below are for the distribution 
P{Sg) of the orbital scattering matrix Sq, in the limit 
of large N. For simplicity of notation we will drop the 
subscript "0" and the superscript "o" . 

For large only moments of the distribution that 
contain as many factors of S as of its complex conju- 
gate S* and in which the indices of factors S and S* 
are pairwise equal are nonzeroiS, This implies that all 
odd moments of P{S) are zero. Since the probability 
of S is close to Gaussian for large N, with small non- 
Gaussian corrections, P{S) is well characterized by its 
moments. The second moment characterizes the Gaus- 
sian part of the distribution; Higher moments represent 
the non-Gaussian corrections to P{S). 

In order to specify P{S), it is sufficient to consider 
generic moments in which all (pairs of) indices are dif- 
ferent. Moments in which two or more pairs of indices 
coincide are fully determined by the invariance of the en- 
semble under basis changes in the leads S For the gen- 
eral time-dependent case, we need the second moment 
{Sij{t, t — T)Sij(s, s — a)* for i j only, which was cal- 
culated in Ref. l48l 



{Sij (t, t - T)Sij (s, s- a)*) 



1 



iVrd 
X exp 



5{T-(7)e{T) 



(B4) 



[ dr' + {a{t -0- a{s - Of + 2 5](X,(i - C) - X,{s - Of^ 



For the time-independent case, we need to go calculate the corrections to the Gaussian distribution. Using the 
short-hand notation S{1) = S'[e(l), Q!(l), Ari(l), ^2(1)], the relevant moments are 



(5., (1)5., (2)*) = P^i(l,2), 

(5,;,(l)5fc,(2)*5H(3)5,i(4)*) = 1^2(1,2,3,4), 

(5,,(l)5,,(2)*5fei(3)5fez(4)*) = M/i(l,2)iyi(3,4)-M^i,i(l,2,3,4), 

(S'y(l)5fe,(2)*5fcz(3)5™/(4)*5™„(5)5,„(6)*) = 1^3(1,2,3,4,5,6), 



(B5) 
(B6) 
(B7) 
(B8) 
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where all indices i, j, k, I, m, and n are different. In addition to the moments listed above, we need moments in 
which the indices of a factor Sij or S*^ are interchanged. These follow from the ab ove moments using the identity 



Sij (e, X_uX2 , a) — Sji {e,Xi, X2 



The function Wi was calculated in Refs. the function W2 was calculated 



in Ref. |48|. Repeating the calculations of Ref. |43for Wi^i and W3, one finds, for large N , 

1 



where 



Wi(l,2) 
W^2(l,2,3,4) 
1^1,1(1,2,3,4) 

^^3(1,2,3,4,5,6) 



^^1(1,2) 



^^2(1, 2, 3, 4) 



Fi(l,2)Fi(3,2)Fi(3,4)i;^i(l,4)' 



Fi(l,2)2Fi(3,4)2 



i^2(l,2,3,4)F2(l,4,3,2) , ^2(1, 2, 1, 4)^2(1, 4, 3, 4) 



4)^^1(3, 2) 



F2(l, 2, 3, 2)^2(3,2,3, 4) ^3(1, 2, 1, 4, 3, 4) ^^3(1, 2, 3, 4, 3, 2) 



Fi (1, 2)Fi (3, 2)i^i(3, 4)Fi (5, 4)^^i(5, 6)Fi(l, 6) 



i^i(3,2) 
^2(1,2,3,4)^2(1,4,5,6) 
J^i(l,4) 



^2(1,2,5,6)^2(3,4,5,2) ^2(3, 4, 5, 6)^2(1,2,3,6) ^ 
H ^TTi^s \ ^ ^^3(1,2,3,4, 5,6) 



F„(l,. 



N 



^(-l)'"£(m)Td 



Fi(3,6) 



^(-l)"a(m) 



i=i 



^(-l)"X,(m) 



(B9) 
(BIO) 

(Bll) 

(B12) 

(B13) 



With the help of these moments, the ensemble averages of Sees. |3 and IVII can be done using the diagrammatic 
technique of Ref. [t^ ^ 



APPENDIX C: QUANTUM DOT WITH 
RELAXATION 

A highly simplified model for relaxation is to consider 
the model of Sec. IIII CI for r < 1. In this appendix, 
we describe a more realistic (but still phcnomenological) 
model for relaxation in the quantum dot. 

In this model, the quantum dot itself (rather than the 
point contact) is coupled to a fictitious lead, see Fig. 
[TUl This fictitious lead is connected to a reservoir, the 
chemical potential of which is chosen such that no net 
current fiows through this lead. For that reason, the fic- 
titious lead is referred to as 'voltage probe'. This model 



was originally proposed by Biittiker,^'^'^^ and has been 
applied by various authors to describe relaxation effects 
in chaotic quantum ciots4^iS4^iS4S& In order to model 
spatially distributed relaxation, the voltage probe has 
M 3> 1 channels, each coupled to the quantum dot via a 
tunnel barrier with transmission probability Fy <C 
The relaxation rate is set by the product MTv- 

Electrons in the physical point contact are labeled "R" 
and "L" as before. We label electrons entering the dot 
from the voltage probe by "V" , and electrons exiting the 
dot towards the voltage probe by "W" , see Fig. The 
scattering matrix of the dot now has dimension (N + M), 
and relates the fields ipu and ipv to the fields ipE and ipw. 



N ^ M 



k=l'' k=l-' 



fc=l'' fc=l' 



(CI) 



The fields ipv and ipw are evaluated at the interface between dot and the voltage probe reservoir. 

In order to derive an effective action for the quantum dot with voltage probe attached to it, we again attach a 
fictitious lead at the location of the point contact, as described in Sec. IIII Cl The Hamiltonian is written as the sum 
H = Hq + Hi , where Hq and Hi represent ideal coupling of the fictitious lead and an impurity at the contact to the 
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FIG. 10: Top panel: schematic drawing ol a quantum dot 
with a voltage probe. Bottom panel: layout of labeling of 
fermion fields "R", "L", "V", and "W". 



fictitious lead, respectively. For Hi we take Eq. (|3.13|) with Vj = 1, j = 1, . . . , N. We then use Eq. (|C1|) to eliminate 
the fields tpL from the expression for Hi and arrive at the effective action 

S ^ 2vF I dti I dr^ [v;^^(0,ii)5^„(ii,ti-T)Vi„^(0,ii-T)+h.c. 
+ 2vF dti j drY^ [V'™l(0, ii)s™p(<i, ii - r)i^^p(ii - t) + h.c. 



(C2) 



7np 



The following formal manipulations closely follow those of Ref. 0. Since the current Ij, j = 1, . . . ,N, does not 
depend on the fields i^vk, k = 1, . . . , M, we can integrate these out. This results in a new effective action 



S' = 2vF I dti I dri^[i/;^^(ii)^,„„(ii,ti-ri)V;„^(ti-ri)+h.. 

+ AVp J^dtidt2 J dTidT2 ^ 1plnLi*l)Smp{h,tl - TijGvih -Ti,t2- T2 ) (s^ )p„ , V'nL (^2) , 



(C3) 



where 



Gvitl - Tl,t2 - T2) = -i{Tc'4>Vp{tl - Tl)'4>\rj,{t2 ~ T2)) 



(C4) 



is the contour-ordered Green functions for fermions coming in from the voltage probe. (Note that the contour ordering 
is with respect to the contour times ti and t2.) We assume that the change of potentials (either external or internal) 
is slow in comparison to the minimum of the relaxation time and the escape rate into the physical leads. Then, using 
the gradient expansion we find 



j dTidT2s{ti,ti - Ti)Gv{tl - Ti,t2 - T2)s''{t2 - T2,<2) = Gy(ti,i2) + 5Gvitl,t2), 



(C5) 



where 6Gv{ti,t2) is defined as 



SGv{tl,t2) ^ ~ j dTidT2S{ti,ti-Ti)Gvitl-Tut2-Ti)S\t2-T2,t2) + ^^ J 



+ 27TVF J de 



X tr 



(dS{t-e) dfiv dS{t;e)\ . (ds[t:e) dfiy ds{t;e)\ . . 



(C6) 
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Here is the distribution function in the voltage probe, and S{t;e) and s{t;e) are the Fourier transforms of the 

scattering matrices S{t,t — T) and s{t,t — T), respectively, see Eq. (I6.1|) . The second term in Eq. HC6(I does not depend 
on the contour positions of ti and ^2- Upon substitution of Eq. (|C5|I into the effective action HC3|I . one finds 



S' = 2vf dti / dri ^ ■0/„^(ti)S'„„(ti,ti - ri)-0„^(ti - Ti) + h.i 



R 



R 



FIG. 11: Left diagram: Schematic drawing of the fictitious 
lead (E) coupling of which to the left-moving fermions through 
the Hamiltonian 1C8II gives the second term of the effective 
action iCTl . Right diagram: Definition of the fermion fields 
R, E, and L after the change of variables of Eq. HC9II . 



At this point, we recognize that the second term in Eq. 
HC7|) is the same one as the effective action that one would 
have obtained by coupling the left-moving fermions to a 
fictitious lead with N channels, with chemical potential 
/iy. This situation is shown schematically in Fig. 1111 
where the fermions in that fictitious lead are labeled with 
the subscript "E" . The coupling that corresponds to the 



second term in Eq. (|C7|) is that of the Hamiltonian 

rn 

+^™sm„L(o)]. (C8) 

In order to eliminate the Hamiltonian IjCSp wc introduce 
new fermion fields according to^ (see Fig. [TT|l 

old new 

i^Emix) i[ll)„iE{-x)9{~x) - llJrnL{x)0{x), (C9) 

where the Heaviside function 6{x) is defined such that 
9(0) = 1/2. For the new fermion fields, the "L" fields 
have the chemical potential fiy of the voltage probe, 
whereas the "E" fields have the chemical potential fj.^ 
of the fictitious lead. Substituting Eq. ljC9p into the re- 
maining terms of Eq. (|C7() , one finds the effective action 



S = VF I dh I [(i^^^(ii)+V>t„^(ti))5„„(ii,ti-ri)V;„^(<i-ri) + h.c. 



(CIO) 



where, as before, all fermion fields are evaluated at the position x = 0+. Although the field ijjEj, j — I, ■ ■ ■ ,N could, 
in principle, be integrated out of the effective action, it is more convenient to keep them explicit and deal with them 
at the level of the perturbation theory in S. No physical observables depend on the distribution function fj^ of the 
field ipE- 

The chemical potential fiy of the voltage probe reservoir must be chosen such as to have no net current flowing 
into that reservoir at any time. For the current ly in the voltage probe we find, in the adiabatic approximation. 



Iv{t) 



ewF 
47r 



deti 



sit; e)Glj,{t; e)s^ (t; e) - s(i; e)G'^y [t; e)S\t; e) - S{t; e)G^j,it; s)s\t; s) 



(Cll) 
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where Dt^^AB = {dA/dt){dB /de) — {dA/de){dB /dt), the superscript "K" refers to the Keldysh component, and 

GB.R{t',s') = -liT^RWU^')), 
GRv{t',s') = -*(Te^^(t')^^(s')>, 

GvR{t',s') = -t{T,i^y{tyj],{s')). (C12) 

Calculation of the Green function GRR{t\ s') can be done with help of the effective action (|(]10p derived above. For 
the Green functions G^vit' , s') and Gynit' , s'), which keep explicit reference to field ipy, a slightly more complicated 
calculation is necessary. After some algebra, we then find 

(Gvfl),,(t,s) = -VFil^jdh y'dTi^Gy(t,ti-Ti),s,L(ti-Ti,ti)(^„^(ti)+7A„^(ti))V;]^(s)y (C13) 



where the average is taken with respect to the effective 
action HCIOI) . 

This concludes the presentation of the formal frame- 
work of the current calculation in the presence of relax- 
ation. The remainder of the calculation proceeds along 
the same lines as in Sec. IIVI The backscattering term in 
the effective action (|(]10ll is equal to the original effective 
action (|3.17(l . but with a sub-unitary scattering matrix 
S instead of a unitary scattering matrix 5*; the remain- 
ing terms in the effective action do not contain the fields 
^R and hence do not contribute to the interaction cor- 
rection to the current. As a consequence, the expression 
for the interaction correction to the current one obtains 
in an expansion in the action (|C10I) is formally equal to 
the expressions obtained in an expansion in the original 
action (|3.17|) . but with a sub-unitary scattering matrix 
S{t, t ~ t) instead of a unitary scattering matrix. 

There is a crucial difference with the calculation of Sec. 
IIVI however: The fact that the action l|C10p features a 
sub- unitary scattering matrix means that, for strong re- 
laxation, the effective action is small, and an expansion 
in the action is justified. This is in contrast to the case 
of a fully coherent quantum dot, where we argued that 
there was no small parameter that justified an expansion 
in the action. We may thus conclude that the interaction 
correction to the dc conductance and the capacitance cal- 
culated in Refs. ^ and'^ describe the case of a quantum 
dot with a relaxation time much smaller than the mean 
dwell time tj. In this sense, although they do not apply 
to the case of a fully coherent quantum dot, Refs. T and 
do capture the first effect of coherent scattering inside 
the dot at temperatures where scattering inside the dot 
is still predominantly incoherent. 

In light of this, the general trend of the results obtained 
in Refs. and is that the charging interaction reduces 
the effective dephasing rate: the weight of scattering pro- 
cesses with a short delay time is increased,^ leading to a 
smaller probability to dephase.^^ (Since the average de- 
lay time must remain constant, interactions must also 



increase the weight of scattering processes with a long 
delay time. If the dephasing time is shorter than the 
average delay time, the net effect of the interactions is 
to suppress dephasing.) A smaller dephasing rate corre- 
sponds to a larger weak localization correction and larger 
conductance fluctuations, as was found in Refs. |3IE A 
detailed analysis of the expression for the current in the 
voltage probe model for the various limiting cases, as we 
did for the fully coherent quantum dot, is of little practi- 
cal use, since the non-interacting (Hartree) contribution 
to the current contains the dephasing rate, which is an 
unknown parameter this model. 
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